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INTRODUCTION 


«  This  report  describes  the  results  of  converting  the  1977  ROBIN*  program 

written  for  the  UNIVAC  1108  computer  (36  bit)  to  a  ROBIN  program  written 
for  the  Interdata  7/32  computer  (32  bit). 

In  making  the  conversion,  it  was  necessary  to  reduce  the  physical 
size  of  the  resulting  object  program  and  limit  the  number  of  input/output 
units  used  by  ROBIN  because  of  the  limited  memory  of  the  Interdata 
computer  (128,000  bytes)  and  limited  input/output  units  (1  magnetic 
tape  unit). 

A  flow  chart  for  each  ROBIN  subroutine  was  made  and  is  included  so 
that  any  unnecessary  calculations  or  programs  could  be  eliminated.  The 
results  of  these  flow  charts  indicate  that  the  subroutine  LEGNDR,  which 
is  a  very  large  program,  (35,000  bytes)  can  be  replaced  by  data  statements 
for  the  cubic  coefficients  PXY1(51),  PXY2(35),  and  PZ2(21).  These  coef¬ 
ficient  values  are  fixed  and  do  not  change  during  the  calculations.  Since 
the  linear  coefficients  change,  during  the  flight,  a  new  program  LINEAR  was 
written  to  compute  the  values  of  PZl  for  points  from  19  to  51. 

All  data  in  the  subroutines  DRAGT  and  ATMOS  have  been  converted 
from  floating  numbers  (full  word)  to  half  word  integers.  This  procedure 
reduces  by  about  half  the  size  of  these  two  programs. 

Comparative  results  between  the  1108  and  Interdata  ROBIN  programs 
are  included  for  benchmark  data  taken  from  an  FPS-16  radar  track  of  a 
.  sphere.  This  data  is  contained  on  a  unilog  tape  library  locator  number 

Y388  at  White  Sands  Missile  Range. 

*For  a  detailed  description  see  references  1  and  2. 
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I.  DESCRIPTION  OF  PROGRAMS 


This  section  discusses  each  subroutine  and  program  used  by  the 
Interdata  system  and  the  UNIVAC  1108  system.  Due  to  the  different 
requirements  of  the  systems  some  routines  are  different.  The  following 
notation  will  be  observed  for  all  routines: 

*  Denotes  calling,  argument  for  the  Interdata  program 
**  This  program  is  not  used  by  the  Interdata  system 
***  This  program  is  not  used  by  the  1108  system. 

A.  ROBIN  Program 

ATMOS  (HI,T,RH00) 

This  program  is  sent  an  altitude,  HI,  (meters)  and 
returns  Temperature,  T,  (degree  Kelvin)  and  density,  RHOO. 
These  returned  data  are  found  by  interpolating  into 
data  tables  corresponding  to  the  1976  standard  atmosphere. 

CONTROL  (GS , N2MID , KSW2 , KSW , IPGE , I Y ) 

Subroutine  used  to  control  the  calculation  of  the 
theoretical  trajectory  and  bias  correction  values. 

CORRS  (N,A,M,B,C2,IS) 

Subroutine  that  calculates  a  correlation  value 
based  on  the  filter  coefficients.  This  constant  is  used 
by  the  subroutine  DEV  to  compute  the  RMS  noise  error. 


DECALT  (JCNT, $ ,  $)  -  (JCNT , IFLAG1)* 

Subroutine  reads  radar  data  until  there  are  5  con¬ 
secutive  decreasing  altitudes  then  selects  the  next  data 
point  corresponding  to  an  even  second  of  time  as  the  initial 
trajectory  point. 

DEV  ( VRHO , VWX , VWY , VP , RHOOO ) 

Subroutine  that  calculates  the  RMS  error  associated 
with  the  winds,  pressures,  and  densities. 

DRAGT  (CD,AMC,RE, IY, IPGE, $)  -  (CD , AMC , RE , IFLAG) * 

This  subroutine  is  sent  a  Mach  number  (AMC) ,  a  Reynolds 
number  (RE) ,  from  these  data  a  three  dimensional  table 
look-up  is  performed  to  find  the  corresponding  coefficient 
of  drag  (CD) . 

DRIVE*** 

Main  program  used  to  read  unilog  tape  and  writes  this 
data  to  unit  14  via  calls  to  RTDATA. 

DRVT  (K,TEMP,RE,AMC,CD,DENTT,HI,LLL) 

Subroutine  calculates  Reynolds  number  and  Mach  number 
based  on  velocity  calculations  and  temperature  and  density 
data  from  subroutine  ATMOS.  The  drag  coefficient  is  then 
obtained  from  subroutine  DRAGT. 

r 

FITON 

Subroutine  called  by  ROBIN  and  TROBIN  to  calculate 
smoothed  velocity  points  by  fitting  a  Legendre  polynomial 
to  the  raw  positional  points.  Acceleration  fits  are  done 
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in  ROBIN  and  TROBIN  themselves.  FIT0N2  is  an  entry  point 
within  FITON  and  is  used  when  the  number  of  points  for  X  and 
Y  smoothing  is  less  than  the  number  of  points  for  Z  smoothing 

INTER  (INTERP , Z2 ,Z3 , IY ,IPGE , STALT) 

Subroutine  that  interpolates  the  one  second  data 
calculated  by  ROBIN  for  altitude  increments  of  either  200 
meters  (INTERP  t  2)  or  1000  meters  (INTERP  =  2). 

LEGNDR  (NPTS ,MPWR,CF,CFSQ,CF1,CF2, NDER) ** 

Subroutine  called  by  SINGLE  to  calculate  Legendre 
polynji  Lai  filter  coefficients  for  degrees  up  to  10  and 
up  to  100  points.  For  both  the  January  78  and  September 
79  ROBIN  programs,  the  cubic  coefficients  PXY1(51),  PXY2(35) 
and  PZ2(21)  are  fixed.  Since  these  values  never  change 
they  are  stored  as  data  statements  in  the  Interdata  program. 
Since  the  linear  coefficients  PZ1  do  change  they  are 
calculated  by  subroutine  LINEAR  in  the  Interdata  program. 

This  approach  saves  32,000  bytes  of  storage. 

LINEAR  (NPTS)*** 

Subroutine  only  used  by  the  Interdata  program  to 
calculate  the  linear  filter  coefficients  PZ1.  See  descrip¬ 
tion  under  LEGNDR. 

READPT (IX,  ICNT , T, SI, S2, S3, $,$)** 

Subroutine  that  returns  a  single  data  point  of  T,X,Y 
and  Z  at  one  tenth  second  intervals.  Program  reads  data 
from  an  array  generated  by  subroutine  READ1. 


READ1  (IX,D,N, JSTAT)** 

Subroutine  Reads  edited  X,Y,Z  and  T  data  from  unit  15 
into  a  24  point  array. 

REAVG  (ICNT,$,$)  -  ( ICNT , IFLAG 3 ) * 

This  subroutine  obtains  5  one  tenth  second  data  points 
from  READPT  and  returns  the  average  to  the  calling  program. 
This  program  in  the  Interdata  system  replaces  subroutine 
READPT  and  RE ADI. 

REAVGT  ($)  -  (IFLAG4) * 

Subroutine  reads  the  theoretical  trajectory  data  from 
unit  3  and  returns  a  five  point  average.  For  the  Interdata 
program  the  theoretical  data  is  obtained  by  a  call  to  sub¬ 
routine  THEOT. 

ROBIN 

Main  program  controls  all  subroutines  to  obtain  total 
ROBIN  output.  Specific  task  is  to  calculate  from  actual 
positional  radar  data,  winds,  temperature,  pressure  and 
density  and  then  add  bias  correction  terms  from  theoretical 
trajectory  to  produce  final  data  output. 

RTDATA** 

Subroutine  reads  one  data  point  from  the  data  buffer 
produced  by  reading  a  unilog  tape.  The  actual  read  is  per¬ 
formed  by  1108  library  subroutine  TAPIO. 

SINGLE  ( IDEG1 .NPTFT1 , CPI , CV1 , CA1) ** 

Subroutine  normalizes  the  filter  coefficients  generated 
by  subroutine  LEGNDR. 
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SLIDE  (IB,ICNT,$ ,$)  -  (IB , ICNT , IFLAGS) * 

Subroutine  shifts  the  data  in  the  raw  data  array 
by  one  second.  That  is  two  data  points  (oldest  time) 
are  shifted  out  the  top  and  two  new  data  points  (most 
recent  time)  are  shifted  into  the  bottom.  New  data  is 
obtained  from  REAVG. 

SLIDET  (IB, $)  -  (IB.IFLAG6)** 

Subroutine  shifts  the  data  in  the  theoretical  raw 
data  array  by  one  second.  Two  data  points  are  shifted  out 
the  top  and  new  data  points  are  shifted  into  the  bottom. 

New  theoretical  data  is  obtained  from  REAVGT. 

TAB  (ESALT , IY) 

Subroutine  prints  a  table  of  amplitude  ratios  for 
both  density  and  wind. 

TABLE** 

Subroutine  reads  from  unit  7  the  biases  generated  by 
the  theoretical  trajectory  and  stores  these  bias  terms 
in  a  bias  array.  In  the  Interdata  program  the  bias  terms 
are  placed  directly  into  the  bias  array  as  they  are 
calculated  by  TROBIN. 

THEOT 

Subroutine  calculates  a  theoretical  trajectory  hased  on 
the  starting  point  apogee  of  position,  velocity  and  accel¬ 
eration  and  falling  of  the  sphere  thru  the  standard  atmosphere 
with  zero  winds.  The  differences  between  the  winds,  temper- 
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atures,  etc.  found  by  TROBIN,  based  on  the  theoretical 
trajectory,  and  data  from  the  standard  atmosphere  tables 
produce  the  bias  correction  terms.  The  program  for  the 
1108  system  writes  the  trajectory  data  to  unit  3.  The 
program  for  the  Interdata  system  returns  the  next  calculated 
point  to  the  calling  program. 

TIFALL  (HI , TFT , 1FT , ZVM, HI2 , COLAPS , IPGE , IY) 

Subroutine  makes  an  inflation  check  of  a  sphere  based 
on  the  time  of  fall  thru  the  various  layers. 

TIFAL2  (HI , TFT , IFT , ZVM , HI2 , COLAPS , IPGE , IY) 

This  subroutine  makes  the  same  check  as  in  TIFALL 
except  for  a  sphere  of  density  .115  or  .165. 

TROBIN  (GS , N2MID , KSW2 , KSW , IPGE , I Y) 

*  Subroutine  uses  theoretical  trajectory  data  to  compute 

density,  temperature,  winds,  etc.  For  the  1108  system, 
trajectory  data  is  'read  from  unit  3  and  the  results  written 
to  unit  7.  For  the  Interdata  system,  the  data  is  obtained 
by  successive  calls  to  THEOT  and  the  resulting  calculations  from 
TROBIN  are  put  directly  into  the  bias  arrays. 

WANGLE  (THETA, WE, WN) 

Subroutine  uses  east  wind  (WE) ,  and  wind  (WN)  components 
and  returns  polar  north  direction  (THETA). 

B.  Data  Editing  Program 

The  subroutines  and  programs  described  here  are  not 
incorporated  into  the  Interdata  program  because  they  would 
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need  to  be  changed  to  accommodate  data  produced  by  sphere 
track  from  a  modified  Nike  Hercules  radar.  Since  the  radar 
has  not  been  modified  to  track  a  sphere,  nothing  has  been 
done  to  the  editing  programs  except  to  better  document 
their  purpose. 

PMR 

Main  program  copies  the  data  tape  onto  a  temporary 
Tile  and  in  the  process  edits  out  some  systematic  dropouts. 
Also  prints  a  list  of  the  data  at  1-second  intervals. 

A 0 

Main  program  calls  the  editing  routines  via  subroutine 
MAIN  and  then  calls  the  main  data  reduction  program  ROBIN. 

MAIN 

Subroutine  reads  the  variables  IDCHECK  and  TSTART, 
writes  them  on  the  output  and  then  calls  subroutine  MDECK 
to  perform  the  editing. 

MDECK 

Subroutine  edits  data  points  by  discarding  those 
points  which  do  not  meet  the  tolerance  requirements.  Data 
thus  rejected  is  replaced  by  interpolating  between  good  data 
points. 

READST  (NT) 

Subroutine  searches  for  the  first  record  label  on  unit 
NT.  Resultant  parameters  for  the  record  are  written  to  unit  6 


READTT  (NT,TI,SI,S2,S3,NA,NCT,NF1,*,*) 

Entry  point  within  READST.  Reads  data  from  unit  14 
until  a  data  point  with  start  time  is  encountered.  Consecutive 
data  points  are  returned  to  calling  program. 

WRTAPE  (NA,HT1,HS1,HS2,HS3,M) 

Subroutine  writes  edited  data  for  input  to  the  ROBIN 
program  to  unit  15.  The  individual  points  are  packed  into 
a  96-word  record  before  being  written  to  unit  15. 


II.  COMPARISON  OF  INTERDATA  AND  UNIVAC  OUTPUTS 


The  wind  outputs  from  the  Interdata  and  UNIVAC  systems  for  the  bench¬ 
mark.  data  is  compared  in  Figure  1.  This  Figure  is  a  graph  of  the  east- 
west  and  north-south  wind  components  at  1-kilometer  intervals  beginning 
at  98  km  and  ending  at  36  km.  Visual  comparison  of  the  data  shows  that 
the  maximum  error  is  1  m/s  in  the  N-S  wind  and  no  difference  for  the 
E-W  wind. 

The  pressure,  temperature,  and  density  outputs  were  not  compared 
graphically  because  er  the  large  dynamic  changes  that  occur  in  these 
variables.  However,  Table  1  contains  the  actual  pressure,  temperature 
and  density  that  was  calculated  by  the  two  systems  at  the  altitudes 
Indicated . 

In  comparing  pressure,  the  maximum  percentage  error  was  0.12% 
which  occurred  at  81  km.  This  percentage  error  is  calculated  by  use  of 
the  following  equation: 


percentage  error  = 


(ABS(PUNIVAC  -  PINTERDATA))  >;  100 
P INTERDATA 


The  nominal  error  for  pressure  data  calculated  for  altitudes  less  than  76  k 
is  less  than  .06%. 

The  maximum  error  in  temperature  is  1  degree.  This  1-degree 
error  occurs  because  only  whole  numbers  are  printed  on  the  output. 
Internally  the  error  would  be  approximately  the  same  as  for  the  pressure 
data.  The  maximum  density  percentage  error  is  0.33%  which  occurred  at 
91  km.  The  nominal  error  for  density  data  for  altitudes  less  than 
83  km  is  less  than  .06%. 

In  all  cases  the  maximum  error  from  Table  1  is  due  to  the  fact 
that  output  printed  did  not  contain  enough  significant  digits.  Using 
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s/w 


from  UNIVAC  and  INTERDATA 


Table  1.  Pressure,  Temperature  and  Density  Comparison  of  UNIVAC  and 
Interdata  Outputs  Utilizing  Benchmark  Input  Data. 


Alt  itude 

(In.) 

Pressure  (MB) 

Temperature 

(Deg  K) 

Density 

(G/M3) 

- aMU - - 

Interdata 

UNIVAC 

Interdata 

UNIVAC 

Interdata 

UNIVAC 

98 

. 00038 

.000)1 

184 

184 

.00072 

00072 

97 

.00046 

.  00046 

176 

176 

. 00091 

00091 

96 

.00056 

.00056 

170 

171 

.00114 

00114 

y5 

.00068 

.00069 

165 

166 

.00144 

00144 

94 

. 00084 

.  00084 

163 

164 

.00179 

00179 

93 

.00103 

.00103 

164 

165 

.00218 

00218 

92 

.00126 

.00126 

170 

171 

.00257 

,00256 

91 

.00152 

.00152 

178 

178 

.00297 

.00296 

90 

.00182 

.00182 

186 

186 

.00340 

,00339 

'•'9 

.00216 

.00216 

192 

193 

.00390 

.00389 

88 

. 00255 

.00255 

197 

198 

. 00450 

.  00449 

87 

.00302 

.00:  02 

201 

201 

.00523 

.00522 

86 

.00356 

.00356 

202 

202 

.00613 

.00613 

t. 1 

.00421 

. 00921 

202 

202 

.00725 

.00724 

84 

.00497 

.00497 

199 

199 

.00867 

.  00866 

33 

. 00588 

. 00588 

198 

198 

.01033 

.01033 

83 

.00697 

.00696 

195 

196 

.01238 

.01238 

81 

.00827 

.00826 

192 

192 

.01496 

.01496 

80 

. 00984 

.00983 

189 

189 

.01806 

.01806 

79 

.01172 

.01171 

188 

188 

.02162 

.02162 

78 

.01396 

.01395 

139 

189 

.02571 

.02571 

77 

.01662 

.01661 

190 

190 

.03039 

.03039 

76 

.01976 

.01976 

192 

192 

.03579 

.03579 

75 

.02)47 

.0234/ 

i  95 

194 

.04193 

.04194 

74 

.02782 

.05781 

198 

198 

.04887 

.04888 

7  J 

.03289 

.03290 

202 

202 

.05662 

.05664 

72 

.05374 

.01874 

206 

205 

.06549 

.06553 

71 

.04549 

.04549 

210 

210 

.07517 

.07524 

70 

.05321 

.05323 

215 

215 

.08583 

.08594 

i  7 

.06199 

.06202 

220 

220 

.09776 

.09889 

68 

.07202 

.07206 

225 

225 

.11130 

.11135 

67 

.08345 

.08350 

229 

229 

.12678 

.12675 

66 

.09629 

.09630 

239 

239 

.14030 

.14014 

65 

.11062 

.11065 

241 

241 

.15944 

.15943 

64 

.12697 

.12702 

242 

242 

.18228 

.18224 

63 

.14573 

.14574 

245 

245 

.20676 

.20672 

62 

.16677 

.16679 

250 

250 

.23183 

.23174 

61 

.19047 

. 19048 

254 

254 

.26049 

.26045 

60 

.21634 

.21675 

254 

255 

.29480 

. 29507 

59 

.24526 

.24510 

259 

259 

.32907 

.32891 

58 

.27898 

.27881 

261 

261 

.37222 

.37204 

Table  1.  Pressure,  Tsaperature  and  Density  Comparison  of  UNIVAC  and 

Interdata  Outputs  Utilizing  Benchmark  Input  Data  (Continued). 


Altitude  Pressure  (MB)  Temperature  (Deg  K)  Density  (G/H3) 


Interdata  UNIVAC 

Interdata 

UNIVAC 

Interdata  UNIVAC 

57 

.31647 

.31628 

269 

269 

.40968 

.40950 

56 

.35816 

.35795 

269 

269 

.46217 

.46191 

55 

.40574 

.40550 

267 

267 

.52765 

.52736 

54 

.46019 

.45992 

263 

263 

.60848 

.60816 

53 

.52233 

.52203 

262 

264 

.68707 

. 68666 

52 

.59285 

.59253 

265 

265 

.77830 

.77794 

51 

.67338 

.67301 

263 

263 

.89124 

. 89088 

50 

.76436 

.76394 

262 

262 

1.01524 

1.01472 

45 

1.44703 

1.44631 

263 

263 

1.91085 

1.90950 

40 

2.78673 

2.78544 

248 

248 

3.90425 

3.90476 

35 

5.52564 

5.52335 

235 

235 

8.17073 

8.16482 

30 

11.40067 

11.39648 

226 

226 

17.53932 

17.51721 

25.140 

28.56268 

28.56044 

140 

140 

70.91200 

70.90540 

data  where  5  (32  bit  floating  point  accuracy)  or  more  digits  are  printed, 
the  error  was  less  than  .06%.  Thus,  the  internal  error  is  less  than  .06% 
in  all  cases. 

The  total  execution  time  for  calculating  and  printing  the  total 
output  was  26  minutes.  Part  of  the  total  output  is  included  from  the 
initial  height  to  an  altitude  of  50  km. 


TES1  OF  ROBIN  BENCHMARK  DATA 
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III.  CONCLUSION 


The  conversion  of  Che  ROBIN  program  from  Che  UNIVAC  1108  syscem  Co 
Che  Incerdaca  7/32  sysCem  has  been  accomplished.  The  resulcanc  ouCpuCs 
agree  Co  wichin  1  m/s  on  winds  and  .06%  on  pressure,  CemperaCure  and 
densicy . 

The  coefficienCs  for  Legendre  polynomial  smoothing  can  be  stored 
in  a  data  table,  thus,  saving  32,000  bytes  of  storage.  However,  in 
storing  the  coefficients  in  a  data  table,  there  must  be  4  corrections 
made  in  the  FORTRAN  source  code. 

The  code  changes  must  be  made  in  ROBIN  and  TROBIN  as  indicated  below. 


ROBIN 
Line  209 
Change  to 
Line  216 
Change  to 

TROBIN 
Line  37 
Change  to 
Line  44 
Change  to 


NEND  =  N2ST 
NEND  =  N2ST 
NEND  =  NXY2 
NEND  =  NXY2 

NEND  =  N2ST 
NEND  =  N2ST 
NEND  =  NXY2 
NEND  =  NXY2 


+  NZ2 
+  NZ2  -  1 
+  N2ST 
+  N2ST  -  1 

+  NZ2 
+  NZ2  -  1 
+  N2ST 
+  N25T  -  1 


Without  these  changes  one  non-existent  data  point  is  smoothed  by 
a  non-existent  coefficient.  In  the  present  UNIVAC  program  the  non¬ 
existent  coefficient  is  set  to  zero,  thus  the  added  smoothed  point 
value  is  zero  and  does  not  change  the  filtered  output.  By  storing  the 
filter  coefficients  in  a  table,  the  non-existent  coefficient  is  not  zero 
and  a  large  error  occurs  in  the  filtered  output. 
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The  same  smoothing  is  done  in  FITON  but  the  correct  code  exists  as 
indicated. 

FITON 

Line  16  NEND  -  NZ1  +  NIST  -  1 

Line  28  NEND  =  NXY1  +  NIST  -  1 

The  program  will  yield  accurate  results  on  the  Interdata  system. 
When  data  from  the  Nike  radar  is  obtained,  further  programs  would  be 
required  to  edit  the  raw  data  before  running  the  program. 
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IV .  FLOWCHARTS 


UNI VAC  1108  GENERAL  DATA  FLOW 


fcvw-- ■  •-  ■  .  or  '•  • 
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UN I VAC  1108  GENERAL  DATA  FLOW 


INTERDATA  7-32  CENERAL  DATA  FI.OU 


INTERDATA  7-32  GENERAL  DATA  FLOW 


READ  IDD  AND  COMPARE 
WITH  ID  IF  THEY  ARE 
t N°T  EQUAL,  STOP.  IF 
YtHEY  ARE,  CONTINUE  / 


SET  IPGE  -  6  AND 
CALCULATE  CONVERSION 
FACTORS 


ADD  CONVERSION 
FACTOR  AND  CORRECTION 
FACTORS 


KLOWCIiAKI  r OR  MAIN  PKOO'KAM  KOiii.'i  -  l.J 


FLOWCHART  FOR  MAIN  F  ROC  RAM  ROBIN  -  1.4 


FLOWCHART  FOR  MAIN  PROGRAM  ROUIN  -  1.5 


flowchart  for  main  program  robin 


IF  HI  LESS  THAU  60  KM 
AND  TIME  AT  MIDPOINT 
IS  GREATER  THAN 
TIME  OF  FALL  TEST  6 
.PLUS  50  SECONDS 


'IF  2  VELOCITY  X 

.midpoint  i  2VC 


RECALCULATE  CONVERSION  FACTORS  ' 
AND  CALL  CORRS  FOR  NEW  CORRECTION 
FACTORS 


f  nzi  =  sf 


CHANGE  LIMITS 
AND  RECALL 
SINGLE 


FLOWCHART  FOR  SUBROUTINE  ATMOS 


FLOWCHART  FOR  SUBROUTINE  BLKDAT 


DATA  TABLES 
FOR  all  common! 
DATA  VARIABLES! 


FLOWCHART  FOR  SLBKOL'TINL  CONTROL  -  1 


43 


FLOWCHART  FOR  SUBROUTINE  CORRS  - 


FLOWCHART  FOR  bUBROL [ I N£  ULV 


46 


flowchart  for  subroutine  dkagt  -  6.1 


4  7 


•  *  .«■“  ft'rt*-  ■  .  *■"  ■*-•***'  ■■• 


FLOWCHART  FOR  SUBROUTINE  LECNDR  -  1 


COMPUTE  COEFFICIENTS 
FOR  ACCELERATION  (SECOND 
DERIVATIVE)  CALCULATIONS 


FUVCHAKT  EUR  SUBROUTINE  REAUFT  -  13.1 


FLOWCHART  FOR  SUBROUTINE  REAVC  -  17.1 


FLOWCHART  FUR  SUBROUTINE  KTDATA  -  14.1 


57 


FLOWCHART  FOR  SUBROUTINE  TAB  -  22.1 


61 


FLOWCHART  FOR  lUBRotriXt  liit-OT  -  24 


f INITIALIZE  VALUES  T0| 
I  ZERO  ANP  SET  TIME 
I  INCREMENTS  TO  .1  SEC 


KLOl.'CiURT  FOR  SU1SKOUTINF  T1FALL  -  23.1 

Tif all  is  dependent  on  I  FT  as  to  where  in  the  subroutine  it  will 
return  to.  The  first  statement  is  a  computed  CO  TO  statement  with 
eleven  possible  cnoices. 

If  IFT  3  1*  return  immediately. 


IFT  »  2 


64 


flowchart  for  SUBROUTINE  TIFALL  -  25. i 


I  FT  =  S 


START 
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FLOWCHART  FOR  SUBROUTINE  TIFALL  -  25.^ 


IFT  -  6 


IFT  =  7,  IFT  =  AND  IFT  *  ARE 
SIMILAR  TO  IFT  *  5. 


IFT  »  9  AND  IFT  -  II  ARE  STRUCTURED 
THE  SAME  AS  IFT  =  3,  BUT  THEY 
HAVE  NUMERICAL  DIFFERENCES 


67 


FLOWCHART  FOR  SUBROUTINE  TROBIN 


FLOWCHART  FOR  MAIN  PROGRAM  PMR  -  El.l 


FLOWCHART  FOR  MAIK  PROGRAM  PMR  -  El.  2 


FLOWCHART  FOR  SUBROUTINE  MDECK  -  E4.5 


SECOND  RETURN 
THIRD  RETURN 


GET  NEXT 
POINT  VIA 
CALL  READTT 


INCREMENT  NUMBER 
OF  RECORDS  READ 


RESET  COUNTER 
FOR  BAD  POINTS 
IN  SEQUENCE 


MARK  15 


WRITE  ON  OUTPUT 
NUMBER  OF  GENERATED 
POINTS 

NUMBER  OF  DISCARDED 
POINTS 

NUMBER  OF  RECORDS 
.  WRITTEN  _ 


INC 

FC 

PC 

:R£MENT  COUNTER 
)R  GENERATED 
)1NTS 

Vwi 

\i 

1TE  ON  15  / 

ENERATED  POINT/ 

INC 

OF 

RJEMENT  NUMBER 
RECORDS  WRITTEN 

/  WILL  THE  X 
'nEXT  GENERATLDN 
POINT  BE  WITHIN 
.001  SECONDS  OF 
THE  PRESENT  > 
X  Data  POINT? / 


INCREMENT 
TIME  BY 
.1  SECOND 


FLOWCHART  FOR  SUBROUTINE  READTT  -  Efa 


FLOWCHART  FOR  SUBROUTINE  READTT  -  E6.2 


84 


V.  PROGRAM  LISTING  FOR  INTERDATA  COMPUTER 


1  0 
.1  1. 
:1 

:1  : 

J.  .4 
:1  • 


.1  8 
i  8 


■:i:i 


$HSSM 

ROBIN 

4  FORT 


C 

C 


PPAUS 

PROG  RIGID  SPHERE  TRRJ 

CROSS 

NO  FlX  3 

NLSTC 

NLIST 


PERL  nr 

DOUBLE  PRECISION  0X1,  DV1,  DZ1,  DRG,  DHMSL,  DHI 


36  NOV  79  R01  S  RFM3 : ROBIN.  FOR 


INTEGERS 
INTFGFR* 
INTEGER* 
INTEGER* 
INTEGER* 
I  NTFGF.R**- 


2  IFLRG,  IFLRG1,  I  El  HU/,  I  FLOGS,  TP,  ON,  I  MX,  JMX,  IF'l 

2  I DEGN1,  I DEGX2,  I  DEG 71 ,  I DEGZ2,  N2H I D 

2  NNV2,  NZ2,  NX1MID,  NX2MID,  N71MID,  NZ2MID,  N2ST 

2  IX,  J,  W,  MIDMRX,  KMRX1,  KMRX2,  NIST,  NXV1,  N21 

2  IBC,  I  BCR,  i  BCM,  IFLRGB,  ID,  IV,  I,  I  DO,  INTERP,  ICNT,  NF'RT 

2  I  CD,  IF'GE,  IBAL,  IFT,  ISM,  I  Ml,  I  HR,  JST  T,  K,  KSM,  KSW2,  NEND 


COMMON  /TE/TFT1,  TFT2,  TFT2,  TFT-1,  TFTG,  TFT6,  TFT?,  TFT8,  TFT 9,  TFT© 
COMMON  XD/X1,  X2,  X3,  VI,  V2,  V3,  71,  72,  Z2X,  23,  Z3X,  C0NX1,  C0NX2,  C0NX3, 

X  C0IIX4,  C0N71,  C0NZ2,  C0NZ3,  C0NZ4,  GS,  IDEGX2,  IDEGX1,  IDEG22,  IDEG71 
X,  N2MID 

COMMON  /TRBL./  I  MX,  EMASHI  <  600!),  B I ASMX  <  60©  > ,  El  I ASWV  <  606 !) ,  8IASR<60O!>, 
X  8 1 RST < 6©G > ,  JMX, B I RSMC < BOO  > 

COMMON  XCONST/RG,  Rl.R,  HMSl. ,  RMS,  ZB,  DIR,  VP,  RM,  GSRG,  RB 
COMMON  /EXTRA/  NXV2, N22, NX1MID,  NX2MID, N71MID, NZ2MID, N2ST 
COMMON  /CORK/  PXV1<51>,  PXV2<35>,  P22<21>,  PZ1<51> 

COMMON  /COR 1 2/  TOSl  ,  TOOL,  TOOL  SZ,  TOCLCZ 
COMMON  /C0N2/  RMK,  ALPHA,  RU,  RRD 

COMMON  /I  NT/  IH,  HI,  BMC,  WF,  WN,  WT,  THETA,  PI  COR,  RHOCOR, 

X  TCOR,  CD,  RE,  VRHO,  VWX,  VWV,  PI  MB,  VP,  VT,  INTBI  T,  ITIM,  CFREQ,  TND1,  TND2 
COMMON  T I ME < 1 60 > ,  X < 1 GO  > ,  V < 1 6© > ,  2 < 1 68 > ,  TM I D < 50  > ,  XM I D < 50 > ,  VM I D < GO > , 
X  2 1*1 1 D < GO > ,  XVM<5©>,  VVM<50>,  2VM<50>,  IX,  J,  N,  MIDMRX,  KMRX1,  KMRX2, 

X  Cl,  NIS  I ,  NX VI ,  NZ1 ,  ZXM  < GO > 

DIMENSION  CP1<100>,  CAICIOOX.  BNRME<2©> 

DRTR  R  NR  ME / 8 OH?- 3 2  TEST  OF  ROBIN  BENCHMARK  DATA 


C 


4© 

41 

42 

•r 

44 
4  G 
48 

47 

48 
48 
GlT 


102  FORMAT <110,  212,  FS.  2,  23 X,  F7.  1,  F6.  1,  F7.  4,  F8.  5,  I3> 

324  FORMAT  < 1 HO,  20R4  > 

325  FORMAT < 68H0  ID  RD  HR  RL R  GS  RG  HMSL  ZB  RMS 

*  DIR  INT/1H  ,  14,  212,  F€.  2,  Fll.  6,  FI©.  1,  F7.  1,  F6.  1,  F?.  4,  F8 
+.  5,  I3> 

329  F  ORMRT < 20R4 > 

400  FORMAT <1H  ,  95X,  27H+*+R.  M.  S.  NOISE  ERROR  IN+*+/3X,  119HTIME  ALT  . 

*  EWIND  NWIND  SPEED  DIR  PRESS  T  DENS I TV  7  VEL  Z  ACC  TDE 

*N  CF  DENS  PRES  TEMP  EW  NW/3X,  120HZULU  METERS  M/S 

*  M/S  M/S  DEG  MB  K  GR/M3  M/S  M/S2  /CC 

/S  **■.♦+*:  P  E  R  C  E  N  T  +■ *  *:  *;  *  M/S  M/S  I1 


86 


51  40:1.  F0RMATC1H1, 95X,  27H***R.  M.  S.  NOISE  ERROR  IN***/3X,  119HTIMF  At.T  . 

52  *  EWIND  NMIND  SPEED  DIR  PRESS  T  DENS I TV  Z  VFL  Z  fiCC  TDt 

53  *N  CF  DENS  PRFS  TFMP  EW  NM/3X,  1 20HZULU  METERS  M/S 

54  *  M/S  M/S  DEG  MB  K  GR/M3  M/S  M/S2  /CC 

55  *  /S  *****RE RCENT*****  M/S  M/S  1 

56  403  FORMAT  <  IX,  12,  2<1H.  ,  I2>,  17,  2F6.  1,  F7.  2,  2<  14,  F9.  5>,  F8.  2,  F6.  2,  2<1X,  E9 

5?  *.  4),  3F6.  1,  2F5.  1) 

58  C 

59  C 

60  IX  5 

61  I  V:r-3 

62  WR I TEC  IV,  479) 

63  479  FORMAT C47H1 JN I VERS I TV  OF  DAVTON  ROBIN  PROGRAM  -  SEPT  197 7  //> 

64  NZ1=19 

65  CALL  L INEARCNZi > 

66  IBC--2 

67  IBCA-2 

68  IBCM-2 

69  IFL  AGB:-  0 

78  ISWr-1 

71  RE=-0. 

72  P I  COR:-0. 

73  RHGMGS-  0. 


74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 


TCOR:-  0. 
TND2=0. 
PI^0. 

ZVC«0. 
COL.APS=0.  8 
NZIMID^IO 
N1ST—17 
NX1M I D=26 
NX2M ID-18 
NZ2MID-11 
KMAX2:-35 
KMAX1 -51 
MIDMAX— 2'6 
N2MID-18 
N2ST— 8 
C2-1. 

NXV1-  51 


91 

NXV2-35 

92 

NZ2-21 

93 

I DE  GX1 —3 

94 

IDEGX2--3 

95 

IDEGZ1-1 

96 

IDEGZ2-3 

97 

1011  FORMAT  <1H8,  12X,  30HXV-VFL 

XV-ACC 

Z-VEl  2- 

■ACC/ 

98 

A11H  -  PTS.  FIT,  418/ 

99 

B11H  DEGREE  FIT,  418 > 

100 

C  THE  FOLLOWING  READS  ARF 

REF'l.ACED 

BV  EQUAL I TV 

STATEMENTS 

87 


.1  fvl 

c 

o 

RFRD  C  7,  779  :■  RNfiMK 

1  0  7 

C 

RERIX7,  1021  ID,  IX,  IV,  ftl  R,  HMSL,  ZB,  RMS,  DIR,  INTERP 

103: 

ID=2777 

1 04 

RLR  32  46 

;1  IV , 

HHSl.  1229 

106 

ZR-  «-V 

10. 

RM9  164  9 

1  08 

DIht  1 

1  Ci4 

•  NTt-KP  1 

1.1  4 

PRO-  8  7.  7977797 

1  1 

RG  f:  71  i.JV  :'.'X  _• 

1  *  *.• 

c 

09  -9  1  0+5.  2885E— 5+  t  S I N < RL  fl/RRD > **-2. > -5.  9E~8+<SIN 

11. 

02.  —9  ':-:l  ,j'- 

1  i  4 

c 

X  RF(  ’>  :*  +  2  .:>  > 

.1  i  l; 

WRI'rf-  1  t.<  774  >RNRMF 

11  6 

WRIT  t  IV,  101.1  > MX VI ,  UXV2,  NZ1 ,  NZP,  IDEGX1,  IDF.GX2,  IDF.621,  I DEG 

11  7 

WRITE <.IV.  7?0>ID,  IX,  IV,  RLR,  GS,  RG,  HMSL,  ZB,  RMS,  DIR,  INTERP 

118 

2001 

FQRI19T 4 1  A  8  ’> 

11  •- 

I CUT  27 

1 

c 

230 

RERDt  1  X  1  ILD 

121 

720 

I D0--7277 

v ?f: 

IF  ■..  IDD-ID  '  7 ''2 •  227,  222 

...  <? . 

333 

JSOT  1 

1 24 

IF':'  TSTT  ED.  ....:•  STOP 

a  3  -* 

GO  TO  226 

1.76 

VB--  <0.  72  677  >  *D  I R +  D I R+D I R 

127 

RB-  C6.  39"  69908  )-+-D  I R+D I R 

178 

ESRI..T'-  6.  6 

1 79 

I  CD-I 

170 

ft  ML  7.  028966 

1  7.:i 

R'.  FOR  -0.  00G1723 

■1 

KU  8.  31437 

■‘i  “* 

RTM 

1 24 

LFl  6.  6, '  '169E+23 

1 3L« 

RSI  OR  8.  3J  .  37 

1 76 

WZ1-  0.  0 

127 

NF'F'T  0 

1  7ft 

PI  MB  8.  0 

1  79 

CZfc  -  COS  1 ZB/RRD> 

140 

’r»  Z  :  3 1 )  <  '■  .7  b»  /  R  H  f 

l'*:1 

TOSL  0  0661 4  78+S l N <  fll  .R/RRD > 

10  7 

7  0"  L  0.  6661 4 78+COS tfll.R/RRD > 

1.42 

T  OCL  Si  -  T OCl  ,+SZB 

1  11 

TOOL  CZ'  TiCl+CZB 

147. 

I  F'ljF-.  X2 

1.46 

DEL-.  27 

1.47 

RM1-  I  JXV1 

1  4ft 

C0NX1'  12.  / <  DEL  +RI41+  <  RN:1  +RN1-1.  >  > 

149 

C 0 N X 2 - C 0 ! ■' X 1 + 7 .  *<2.  ♦RNl+RHd  -7.  > **?/<  DEL*AN1* < AN1+RN1-1.  > 

1.00 

d  *tftNi*ftX1  -4.  > * < RN1+RN1  —9.  >  > 

88 


1  51  RN1=N21 

152  C0N21-12  /  <  DEL  *RN1*  <  RN1*RN1— 1.  >> 

153  C0NZ2  C0N21+7.  ♦•<3.  *RN1*RN! -7.  >**2/<DEL*RWl*<RNl*RNl-l.  > 

154  l*<RNl*RNl-4.  )*<RNl*flN1-9.  >> 

155  RN1- -NXV2 

3.56  DEL'  1. 

157  C0NX3— 12.  /  <  DEL+-RM1*  <  RN1*RW1— 1 .  >> 


158 

C0NX4  C0UX3+7.  *<3.  *RNl*RNl-7.  )**2/<DEL*RNl*<RNl*RNl- 

159 

l*<RNl*RWl-4.  >*<RNl*RNl-9.  >> 

160 

RW1---N28 

16:1 

C0NZ3-12.  /  <  DEt  *RN1 *<RN1*RN1-1.  >> 

162 

C0NF4-C0N23+7.  *<3.  *RNl*RNl-7.  >**2/<DEL*RNl*<RNl*RNl- 

163 

l*<RNl*RWl-4.  >  *  < RN1+RN1-9.  >> 

164 

CRI  .L  CORRSCNXVl,  PXV1,  NXV2,  PXV2,  RHSM,  2) 

165 

CR!  L  CORPS <N23 ,  P21,  NZ2,  P?2,  RMSD,  2) 

166 

C0NX3-C0NX3+RMSW 

167 

C0NX4= C0NX4+RMSW 

168 

C0N32-  C0N23+RMSD 

169 

C0N24-C0N24+RMSD 

170 

VP-0.  0 

171 

TFT  -0. 

178 

TFT1M3. 

173 

TFT2-0. 

174 

TF  T3=0. 

175 

TFT4-  0. 

176 

TFT5“0. 

177 

TFT6:-0.  0 

178 

TFT7=0.  0 

179 

TKT8-0.  0 

1 80 

TFT9::  0.  0 

1 81 

TFT0—0.  0 

182 

IFCRBSCDIR-1.  X  LT.  0.  02.  RND.  BBS  CRMS-.  115).  LT.  0.  01  )  GO 

183 

IBflL-i 

184 

IF<flBS<DIfl-l.  X  LT.  0.  02.  RND.  RBS<RI1S-.  165).  LT.  0.  01  >  GO 

185 

WRI TE< IV, 410) 

186 

410 

F ORMRT < 20X,.  46HBRI  L  OON  NOT  ROBIN.  TIME  OF  FBI  L  TEST 

187 

IPGF:-IF'OE+l 

188 

IFT-1 

189 

GO  TO  230 

190 

299 

IBRL-2 

391 

306 

I  FT- '2 

198 

230 

CRL.L  DECRL  T <  ICNT,  IFl  RG1> 

193 

I F  <  I FL  RG1 )  301 , 35,  90 

194 

301 

WRITbXIV,  4O0) 

195 

IPGE-IPGF+3 

196 

INTRLT-0 

197 

1001 

DO  52  J-l, KMRX1 

198 

52 

CRI.L  RERVG<  ICNT,  IFLRG3) 

199 

I F  <  I FL  RG3  )  1 003,  40,  90 

200 

1003 

I F  <  KMRX1.  NF.  NXV1 )  GO  TO  3 

TO  300 


TO  299 


USED 


r 


* 


TOO 


20F. 


I 


’  •  t 


KSW2=t> 

GO  TO  t 
:<  KSN?-  2 

!  1 F  (.  KMRXV  FQ.  NXVi'i  GO  Til  3 
Foil  2 
GO  TO  6 
0  K"0  3 

6  1 F  <  T  FLRGB  FO  ©>  CRN.  COMTRI  <GS,  N2T1ID,  KSM2,  KSI4,  IPGE,  IV> 

IFLRGB--  1 
GOTO*,?..  91,  KSt-12 
3  L  •  S  M  KMRX2 
0  K 

C2U.  F I  "ON 

CRL<_  Si.  IDfr  <2,  I CNT ,  I FLAGS > 

IFcilFL  fluG ;■  8,  46-  96 
F  CONTI HUE 
CO  TO  1 

9  L'O  -i'.1  K  i...  t-  *■  sRX2 
N  -  K 

CRi  !  F  7  TO  HP 

.  r  :Ll  Si .  i  OF  C  2  ■  I  CNT,  I FLAGS  > 
l  F  1  l  r  I  ALT  >  ■  1  O,  TO,  OO 
l_i  0:  .i  “  I  Hi  _ 


V  C.  L 

Zj  iO.  0 

GO  TO  .:.j  il.  o  KSM 


2_  t  ■'  i 

4  ;2 

Lo  j  3  £  •  •  1 . HXV2 

: :  2::>pxv2'  I  •*zxii< I 

'  '  -‘-F-:  '-'j'  I  +XVM  <;  I  :> 

•.:!  _:•  >  ' 

i 

....  -■•MVOT  I  >*VVM<  I  > 

:  ’■  NO  H2ST+NZ2  ~± 

?  >. 

00  14  I  M2  ST,  NF.MD 

•  “•  r  • 

r  I-N*’ST-»  x 

IT 

Z2  -  ?.  ?  +P"  :  C  K  #Z VM  <  I  > 

.  •' 

GO  "C  12 

t  !-! 

13 

C  O  22  1,  M72 

‘I 

Z3-2:-.-*.-Z2<  I  >*ZVM<I> 

24iri 

2UT2  rMV2+N2SY  -1. 

24. l 

CO  J.r  I  N2ST,  MFMD 

*2  '"K  ' 

Y>-  I-N281  +1 

,"'4  3 

1201  ZOK+F'NVi!  K  > *2XM F  I 

244 

X3=  :  «  t-.-2.V2  < K > *XVI1  <  I  > 

2 'i3 

-i  "7 
■A.  1 

V3-:  V 2+p:-:V2 (  K  > +WM <  I  > 

18 

:  ;3--X3t'C2 

,:1? 

V Jr  V3+C2 

24 

Z'-.'A  ?3X*C2 

24'j 

?  j.  ■'_:*C2 

.’Or, 

:>■■'•■  VVH<N2|V|ID> 

90 


V2=VVM<N2tlID? 

Z2X=  7XIKN2MID? 

Z2ZV‘I1<N2MID? 

XI -XM I D  <  N2M ID? 

V1--VM I D  <  N2M I D  > 

21= ZMID<N?MID? 

TIM  TM I D  <  N2M ID) 

HI==Zi+<Xl*Xl+Vl*Vl>/<2.  *RG?+HMSL. 

IH  HI 

FNDH  I  ZM I D  < 1  >  +  <  XM I D  <  1  )  **2+VM  I D  <  1  ?  **2  ? / <  2.  *RG  >  +HMSL 

OF'HGR:  1  d+HI/RG 

MTHR-GS/  <  RG*OPHOR*OPHOR*OPHOR  ? 

IF  <  TF  Tt.  EG.  6.  9.  HUD.  HI.  LT.  60000.  ?  TFT6--TIM 
IF < TFT  1  NF.  6.  .  OR  HI.  GT.  80606.  >  GO  TO  500 
TFT1--TIM 

506  IF TFT2.  NF..  0.  .  OR.  HI.  GT.  70660.  ?  GO  TO  561 
Tl-  T2-TIM 
TKT:  TFT2-TFT1 

501  IFCTFT3.  NF.  0.  .  OR.  HI.  GT.  60660.  ?  GO  TO  562 
TFT3  Till 

TFT :  TFT3-TFT2 

502  IFCTFT4.  NF  6.  .  OR  HI.  GT.  55000.  >  GO  TO  507 
TFT4--T  JM 

TFT=  TFT4-TFT3 

503  IFCTFT5.  NF  0.  6.  OR.  HI.  GT.  50600.  )  GO  TO  505 
TFT5-  TIM 

TF  T :  TFT5-TFT4 

565  IFCTFT7.  NE  0.  6.  OR.  HI.  GT.  45606  >  GO  TO  566 
TFT7=  TIM 

TFT  TFT7-TFC5 

566  I F < TFT 8.  NF.  0.  0.  OR.  HI  G1 .  40606.  ?  GO  TO  567 
TFT8  Till 

TF  T -  TFT8-TF  T7 

567  IF <TFT9.  NF.  6.  6.  OR.  HI.  GT.  35660.  ?  GO  TO  568 
TFT9-TIM 

TFT--TFT9-TFT8 

508  IFCTFT0.  NE.  0.  0.  OR.  HI.  GT.  30060.  >  GO  TO  564 
TFT0-TIM 
TFT:  TFT  6-TF1  9 

504  IFc  IBHl  .  EG  1?  GO  TO  569 

Cftl.l.  TIFRI  .2  <  FNDH  I,  TFT,  IFT ,  ZVM,  FSRLT,  COLRPS,  IPGE,  IV? 

GO  TO  51 6 

5G9  CHI  l  TIFHI  l  <  FNDH  I,  TFT,  IFT,  ZVM,  ESRt  T,  COLRPS,  IF’GF,  IV? 

5.16  TFl  -0. 

IF  CHI.  GT  BJRSHICl??  GO  TO  39 

FNDLO-  Zll  I D  <  KMRX2  ?  +  <  XM  I D  <  KMRX2  >  **2+YM  I D  C  KMRX2  ?  **2  ?  /  < 2.  *RG  ?  +HMSL 
C  0  R  X :  -  T  0  C  L  S  Z  +■■  Z  2  X + T  0  S  L  V  2 
CORV::-  T 0CLC2*Z2X— TOSL  ♦XZ' 

COR?  — T0CLSZ*X2-T0CLC2*V2 
GO  TO  <19,  38?,  ISM 


91 


Hi  •/•HI.  T,  RHOU  3 

Pi.*  7_7'  +  !i‘  I  If-T  1  1  0+  XX+X_3  +  V_i.'+-V.£ 3  /  K  ZX- WTHR*  C  71  +RG 3  3  :+;+2  3 
I  hi  ;  h'Hui.  i+VT  :+  '  T +110.  4  }  /  (  0.  00H001  4  5 8 -3 '"14’ T  c!  T + T  -+: T  3  3 
’  >:.  :',<R7  2&Q.  16XT3  3/348.  29205 
!  ‘rtfj  i  I  1  CO-  Hl'iC,  RE.  IF(  AG  3 
i.  rlG  'Gil,  2000,  51  1 


i .  !  CT  i ,'.060'  GO  TO  20 

Hi  •  ?V! ! J  > /?  VM •' NZ2 3  >  > / <  ZM I D  < 1 3 -ZM I D < NZ2 >  3 *2  © 

7  !  b  ■  - :  ■ :  i  I  t  .  ,J  :  3  I.' 

.0bi-4— bL'R  f  .  _-'5E— 8  +  G07197EI2V  DEI  23 
0O03  *  +  2.r'C-T  .•  -•rir_:;;+  0Oi7j  9732, •••'DFI  .2.J 

:  ,h  L  f  i  v  13  :•  GO  TO  47' 

'  77  :■  :i  L  r.  R."  1  3  GO  TO  20 

-  i  [  V  ■  •  '  *7  0  3 

•■.  1  GR-IDENSnv  GRADIENT  EXCEEDS  NOMINAL  VALUE  > 

„+■'  c+i 

IP'S  ~0  0.  3  COLAPS'- 1 H 
.  ,  Q  -yr;  4.RH0 /AMS  >  +WTHR 
X  Xi+WTMO 
r  V;1  31.  IT  i  ID 

,  V  .  rj..,  'l  'Q 

i  _  - '  7  iFt'  < — BuUVZ 

"7-i,r  -  ;•■•■:  7 + C  0  R  X  -  P  0 1 J  V  X  >  W  DE  N  0  T1 

-  ;  < vs+corv-bguv v > /wofnom 

.  f-fI+RG’3 

-IS'-*  '  1  -.4.7 .3  OF.  0.  1,3  GO  TO  22 

-  i 

.  ;  '7-0X3  3t2+  <  V 7 - W V 3  +  + 2  +  <  7  2 — W Z 3  +■ * 2 3 

f'i ’  3  '  O0E7~  7  .i  +  |,.jTHR+  i.  Rg+71  3  3  /  ‘  I.X3-4'AB;+;lv,T;+-  C  Z2— MZ3  +WTHR  +  *•’.  RG+  Z1  3  +  V.'E'  3 

-I',.  L.E  O  3  GO  TO  80 

:•  .  : ,  24  3 ,  r  sm 

3PO+PHU  .'H  M  T 

.  j.  *-G'-  .3'r  .G  t  t  —  >  .  H I  —  H 1 1  3  +  EXP  (G.  5*01  _0G < RHO+RHOOQ  3  3  /  <  T  RG+0.  5*<HI+HI  I 
’  RU+  RHiJ  3 

.  3  SnpT  ''288  :1*/T3  3/340  29265 

lio  i'7  3  'v'  i  3  c  T  +II0  4  3 /  v  0.  0G6061458*SQRT  7  T '3  3  3 

1  Phut  COO,  A  TIC..  RF,  I  FI  AG  3 


'•rv  ■'  •+.  mmy 


IF  <  I  FI.  RG>512,  2m&,  51  2 
512  I  CD  1 

IFCMPRT.  EQ.  0)  GO  TO  26 
CO~ CDO+  <  CD-CD8  ) X3. 

I F  flF:S < RHO-RHOO ) /RHO.  LT.  ©.  003333)  GO  TO  28 
I FL  NF'RT  GE.  5>  GO  TO  27 

26  RHOO-RHO 
NF'RT  '  NF'RT+1 
GO  TO  2G 

27  RHO:  C  RHO+RHOO  )  /2.  0 
C  D  -  <  C  D + C  D  0  )  /  2 .  Q 
NF'RT-  0 

GO  T  0  2G 

28  RHO 00- RHO 
NF'RT :  0 
RHOO  RHO 
TBS:  0.  0 

R  Ei  — . ;  -  0 .  01 
WXBS  0.  8 
WYE'S  0.  0 

I F  I BC+1  GE  JMX)  GO  TO  45 
IF'.: HI  I  T  60000.  >  GO  TO  45 

IFCRI'iC  LE.  RI1C1.  AND  AMC.  GE  .  85.  AND.  RHC.  LF.  1.  3)  GO  TO  460 

46  IFCBIASHICIBCR).  LE.  HI)  GO  TO  47 
IBCA:-  IBOR+i 

GOTO  46 

47  I HI —I BOA 

I BO- I BCR-1 

RAT  10'  <  H I  -  B I  ASH  I  <  I  111 )  )/<BIASHI  <  IBO-BIRSHI  <  IM1)  ) 

GO  TO  475 

4  60  I F  <  B I RSHC  < I BCM > .  L F.  AMC )  GO  TO  470 
IBCM'  IBCM+1 
GO  TO  460 
470  I  HI:- IBCI'1 
I  BO:  I BCH—1 

RRT I  O'  <  AHC-B I ASMC  <  I  111 )  )  /  <  B I ASHC  <  I BC  )  -B I ASMC  <  I  Ml )  ) 

475  TBS-  BIRSTC  I  I'll ) -K  B I  AST  <  IBO— BIAS  f<  1111)  > ♦RAT 10 
RBS--B I RSR  <  I  111  >  +  (BI  RSR  <  I  BC  )  -B I ASR  <  I  Ml )  )  ♦RRT  1 0 
FIXES'  B I ASWX  <  I  111 )  +  C  B I ASMX  <  I  BC  )  -B I ASWX  <  I  Ml  >  )  ♦RAT  1 0 
HVBS'  BIR'SWVC  IMD  +  CBIASWVC  IBO-BIASWVC  III  t  )  ) ♦RAT  10 
AI1C1  -RI1C 
45  T COR  T-TBS 

R  HOC  OF'  RHU— RBS 
W  X  C  U  F' ;  W  X — W  X  B  -• 

WVCOR  WV-WVBS 
P I  COR:  T  CM  R  ♦R  U  ♦  R  H  0  C  0  R  X  A 1 1 K 
W  E = W  X  C  0  R  ♦  S  7  B + W  V  C  O  R  *  C  Z  B 
UN-  W  X  C  O  R  ♦  C  7  A  -  W  V  C  O  R  ♦  S  2  B 
WT  '  SORT  <  WF**'2+WN*+=2) 

CALL  HANOI  CTHFTA,  WF,  WN)  Q, 


401 

IHETA-THLTA 

462 

RHOMGSr  RHOCOR*1000. 

46.? 

CALL  DEV  <  VRHO, VWX,  VWV,  VPR,  RHOMQS  > 

464 

VP=SQRT  <  VF’R  >  /  <  P I  *999.  831  > 

405 

VT-SGRT  <  VP**2+VRH0+*2  >  *100.  0 

406 

VP=VP*100.  0 

4H7 

VRHO-  VRHO>1©0. 

408 

Vi-4  SORT  C  VWX+C2F5  >  **2+  <  VWV+S2P  >**2> 

469 

VE=  -SORT < C  VWV*C2B > **2+ <  VWX+S2B >  **2 > 

410 

VWX- VE 

411 

VWV- VN 

41.7 

TND2--NA+P  I COR 7  < RSTAR+TCOR > *1.  0E-6 

41"? 

1  ri..  -TIM 

414 

H I Nr  1 TIM/60 

415 

I  SEC-- 1  TIM  -60+MIN 

416 

I  HR'  1-1  IN/60 

417 

MIN^-MIN-60+IHR  . 

41B 

PIM;-l-  0.  61+PICOR 

41.9 

CFREG:  PICOR+7.  80064E+5 

420 

TNDl-  RHOMGS+RTN 

42:1 

IT-TCOR+.  5 

4  c-.' 

I R  Ir.  •■-RE 

423 

PIT  PI 

424 

HI  I --HI 

425 

ALT;  IH 

426 

I F<  INTERN.  GE.  1>  CALL  INTFRUMTFRP,  22,  23,  IV, 

IPGF,  B I  ASH  I  <1  >  > 

427 

IFCIPGE.  L.T.  56 >  GO  TO  33 

428 

IPGF-- 3 

479 

hr: tec: iv,  401  > 

430 

-*:•  •• 

CONTINUE 

431 

TFCHJ.  GT.  BIASHI<1>>  GO  TO  43 

432 

IF <  INTFRP.  EG.  0>  WRITE <4,  3000 >  ALT,  T,  RHOMGS, 

PIMB,  WE,  WN 

3000 

FORMAT  <  '  6E10.  4> 

434 

ISM- 2 

435 

I  Fc;  I  NT  ERF.  FQ.  2>  GO  TO  2060 

436 

IPGF- -  IPGF+1 

437 

WRI  iE(IV.  403 >  Ii  IR,  MIN,  ISFC,  IH,  WE,  WN,  WT,  IHETA,  PIMB,  IT,  RHOMGS. 

438 

+TMD2,  CFRF.Q,  VP. HO,  VP,  VT,  VWX,  VWV 

429 

GO  TO  2006 

440 

43 

CONTINUE 

441 

IFCIMTERP.  EQ.  2)  GO  TO  2050 

442 

IF'GF:  -IPGE+i 

443 

WR I TE <  IV,  434 >  I HR,  MIN,  I SEC,  IH,  22,  23 

444 

434 

FORMAT (.  1 X,  12,  2<1H.  ,  I2>,  17,  45X,  F8.  2,  F6.  2> 

445 

GO  TO  2050 

446 

2000 

ICD-ICD+1 

447 

IF<  I CD.  LE.  5>  GO  TO  2050 

448 

I  CD-I 

449 

ISM-1 

450 

2050 

CONTINUE 

94 


451 

NPRT  *0 

452 

32 

K=KMAX2-1 

453 

DO  31  1-1/  K 

454 

IP1-I+1 

455 

XMID<I)=XMID<IP1> 

45ft 

VMID<I)=VMID<IP1> 

45? 

2MID< I >=ZMID< IPl) 

458 

TMID<  I  >~-TMID<  IPl) 

459 

XVM< I >-XVM< IP1> 

460 

VVMC I )=WM< IPl) 

461 

ZXM<I)=ZXM<IP1) 

462 

31 

2VM< I >=ZVM< IP1> 

463 

GO  TO  <23/  30)/  KSN2 

464 

29 

N-KMRX2 

465 

CALL.  FI  TON 

466 

CALL.  SLIDE<2/  ICNT/  IFL.RG5) 

46? 

I F  <  I  FLAGS  )  240/  40/  98 

468 

3© 

N-KMRX2 

469 

CALL  FIT0N2 

4?0 

CAl  L  SLIDE <2/  ICNT/  IFL.AG5) 

471 

IF<  IFL.AG5)240/  40/  98 

472 

240 

IF<HI.  L.T.  60800.  .  AND.  TIM.  GT.  TFT6+50.  )  GO  TO  250 

473 

GO  TO  260 

474 

250 

IF  <NZ1.  EG.  51  )  GO  TO  260 

475 

NZ1-NZ1+4 

476 

TFT  6-TFT  6+50. 

47? 

NZ1MID=NZ1X2+1 

478 

N1ST:=  I ABS  <  NX1M I D-N21M I D  )  +1 

479 

2  VC-2VI1  <  KMAX2  ) 

480 

CAl  l.  LINEAR<N21) 

481 

GO  TO  11 

482 

268 

IF  <2VM<N2MID).  NE.  2VC)  GO  TO  11 

483 

AN1-NZ1 

484 

DEL-.  25 

485 

C0N21  =12.  /<DEL.*AH1*<AN1*AN1-1.  )) 

486 

C0N22-C0NZ1+7.  *<3.  *ANl*RNl-7.  )  *+2/  <  DEL*AN1*  <  AN1  *AN1 

487 

1*<AN1*AN1— 4.  >*<ANl*ANl-9.  )) 

488 

CQNZ3-C0NZ3-RMSD 

489 

C0N24=C0N24-RMSD 

498 

CALL  C0RRSCN21/  P21/  N22/  PZ2/  RI1SD/  2) 

491 

C0N23=C0N23+RMSD 

492 

C0NZ4-C0NZ4+RMSD 

493 

GO  TO  11 

494 

C  40 

READ  < I X )  IDD 

495 

40 

IF<  ID-IDD)  90/36/98 

436 

36 

ICNT-25 

497 

INTAl.T-0 

498 

IF<HI.  LT.  50000.  )  GO  TO  381 

499 

IF<HI.  LT.  60800.  )  GO  TO  1501 

508 

IF<HI.  LT.  70000.  )  GO  TO  1502 

95 


501 

TF<H1.  I  T.  8S000.  >  00  TO  1503 

502 

GO  TO  301 

503 

1501 

IF1  -1 

504 

GO  TO  301 

505 

1502 

IFT-10 

506 

TFT3—999. 

50? 

GO  TO  303 

50:-: 

150": 

I  FT  8 

509 

TFT2-  999. 

51  0 

GO  TO  SOI 

51 1 

C  25 

READ  v I X >  IOD 

51 

~>r— 

I F  C  •;  D-  :  CO  ;•  1 500..  37,  1 56G 

5  ..*•• 

ICN‘>  25 

514 

GO  TO  220 

51 '  ■ 

©o 

WRITE* IV,  89  > 

53  b 

89 

FORMAT <  '  i# *■*■*■*  DENSITY  IS  NON-POSITIVE 

51? 

90 

h!  T =-999. 

51 

hr 1 1  r  :  tv,  93  > 

53  9 

91 

mjKMAT  ■  /2/10X,  19HN0RMAL  END  OF  ROBIN> 

520 

WRITfc  <4,  3O00 >  Rl  T,  T,  RHOMGS,  PIMB,  WE,  WN 

521 

WRITE  <4, 3O00>  COLRPS 

tr 

OV.  l  tab<esri  t,  iv> 

r»r*2 

1  TUT  -999 

-i*.  ^  i  i  1 1  i 

524  1 500  STOP 

4RSSI-1 


EXECUTION  HRLTED  *****'  ) 


4-FORT 


l  1ST 


$HSSM 

ATMOS 

*FORT 


SCRAT 

PROG  TF.MP&DENSE  CALC 

CROSS 

N0RX3 

NLSTC 

NL.IST 


28  NOV  ?S>  R01  S  REMS :  RTMOS.  FOR 


SUBROUTINE  RTMOSCHI,  T,  RHOO> 

INTE0FR*2  TP <302},  DN<307},  J,  TNPW,  I..  KH,  J0FS<5> 

RERL  HSTS<5},  HSTP<5} 

RF VISED  TABLE  FOR  1976  STANDARD  ATMOSPHERE 
50KM  -  1GOKM  IN  8.  5KM  STEPS 
100KM  -  141 KM  IN  1.  ©KM  STEPS 

DATA  JOFS,  HSTP/261,  161,  71,  51,  1,  1000.  ,  500.  ,  280.  ,  180.  ,  208.  / 
DATA  HSTS/10000'0.  ,  50000.  ,  32000.  ,  30000.  ,  20000.  X 
DATA  <TP<I>,  1-1,  74}/  21665,  21678, 


28 

X 

21698, 

21718, 

21738, 

21758, 

21778, 

23797, 

21 

X 

21817, 

23  837, 

21857, 

21877, 

21897, 

21917, 

22 

X 

21936, 

21956, 

21976, 

21996, 

22016, 

22036, 

23 

X 

22056, 

22075, 

22095, 

22115, 

22135, 

22155, 

24 

X 

22175, 

22194, 

22214, 

22234, 

22254, 

22274, 

25 

X 

22294, 

2231 3, 

22333, 

22353, 

22373, 

22393, 

26 

X 

22413, 

22432, 

22452, 

22472, 

22492, 

2753  2, 

2? 

X 

22532, 

22551, 

22571, 

22593 , 

2261.1, 

22633, 

28 

X 

22658, 

22660, 

22678, 

27680, 

22690, 

22700, 

29 

X 

22710, 

22720, 

22730, 

22740, 

22750, 

22759, 

38 

X 

22769, 

22779, 

22789, 

22799, 

22809, 

22819, 

31. 

X 

22829, 

22839, 

22849, 

22875, 

22931, 

22986/ 

32 

DATA  <TP<  I  >,I-* 

75,  160}/ 

33 

X 

23041, 

23097, 

23152, 

23208, 

23263, 

23318, 

34 

X 

23374, 

23429, 

23485, 

23540, 

23595, 

23651, 

35 

X 

23706, 

23762, 

23837, 

23872, 

23928, 

23983, 

36 

X 

24039, 

24094, 

24149, 

24205, 

24260, 

24315, 

37 

X 

24371, 

24426, 

24481, 

24537, 

24592, 

24647, 

38 

X 

24783, 

24758, 

24813, 

24869, 

24924, 

24979, 

39 

X 

25835, 

25090, 

25145, 

25200, 

25256, 

2533  3 , 

48 

X 

25366, 

25422, 

25477, 

25532, 

25587, 

25643, 

41 

X 

25698, 

25753, 

25808, 

25864, 

25939, 

25974, 

42 

X 

26029, 

26085, 

26140, 

26195, 

26250, 

26306, 

43 

X 

26361, 

26436, 

26471, 

26526, 

26582, 

26637, 

44 

X 

26692, 

26747, 

26802, 

26858, 

26913, 

26968, 

45 

X 

27023, 

27065, 

27065, 

27065, 

27065, 

27065, 

46 

X 

27065, 

27065, 

27065, 

27065, 

27065, 

27065, 

47 

X 

27065, 

27065/ 

48 

DATA  <TP< I },  1^ 

163, 229}/  27065  ,  27065  ,  2706! 

49 

X 

26903, 

26765, 

26628, 

26498, 

26352, 

26215, 

50 

X 

26077, 

25940, 

25802, 

25664, 

25527, 

25389, 

97 


5:1 

X 

25252, 

25:114, 

24977, 

24840, 

24702, 

24565, 

52 

X 

24425, 

24290, 

24153, 

24015, 

23878, 

23741, 

5:? 

X 

23604, 

23466, 

23329, 

23192, 

23055, 

22918, 

54 

X 

22781, 

22644, 

22506, 

22369, 

22232, 

22095, 

55 

X 

21958, 

21321, 

21685, 

21548, 

21426, 

21328, 

56 

v 

n 

21231, 

21133, 

21035, 

20938, 

20840, 

20742, 

5? 

X 

2Q645, 

20547, 

20449, 

26352, 

20254, 

20156, 

53 

X 

20059, 

19961, 

19864, 

19766, 

19669, 

19571, 

5.9 

X 

j  9474, 

19376, 

19279, 

19181, 

19084 

7 

r.0  DFlTft  <TPCI>/  I  "230,  302  >7  1393?, 


/.) 

V 

18339, 

18792, 

18687, 

18687, 

18687, 

13687, 

tv;' 

y 

18687, 

18687, 

18687, 

18687, 

18687, 

18687, 

18687, 

18689, 

18696, 

18708, 

13725, 

18747, 

6*1 

»•*» 

18774, 

18805, 

18842, 

18384, 

18931, 

18383, 

65 

X 

19040, 

19104, 

19172, 

19247, 

19328, 

19415, 

66 

rt 

19508, 

19716, 

19953, 

20223, 

20531, 

20884, 

6 

X* 

21.289, 

21763, 

22329, 

23033, 

24000, 

25200, 

66 

7, 

26400, 

2760©, 

28800, 

30000, 

31200, 

3240  , 

6S* 

u 

ft 

3360, 

3488, 

3600, 

3719, 

3835, 

3949, 

70 

V 

4862, 

4172, 

4280, 

4386, 

4430, 

4592, 

7:1 

x 

4693, 

4791, 

438?, 

4982, 

5075, 

5166, 

-I.- , 

«*  »** 

ft 

5255, 

5343, 

5429, 

5513, 

5536, 

3678  7 

C 


OfiTR  <DMC  I  >,  I~-l,  74  >7  8891,8611, 


r 

v* 

8338, 

8074, 

7818, 

7571, 

7332, 

7181, 

i  r< 

I*s 

6877, 

6660, 

6451, 

6248, 

6052, 

5862, 

r*  i* 

X 

5678, 

5580, 

5328, 

5161, 

5000, 

4844, 

»■  o 

X 

4693, 

454  7, 

4486, 

4269, 

4126, 

4008, 

79 

X 

3884  , 

£  i 

3647, 

3534, 

2425, 

3326, 

80 

X 

3217 , 

3118, 

3022, 

2929, 

2839, 

2752, 

8:1 

X 

2663, 

2586, 

2507, 

243:1, 

2356, 

2284, 

'' 

X 

2215, 

2147, 

2882, 

2©1 9, 

195 7, 

1898, 

x 

1841, 

1812, 

1785, 

1758, 

1731, 

1704, 

84 

X 

1679, 

1653, 

1628, 

1603, 

1579, 

1 555, 

85 

X 

1531, 

1.508, 

1485, 

1462, 

1440, 

141  9, 

86 

X 

1397, 

1376, 

1355, 

131 4, 

1273, 

12337 

84 

DftT H  < 

DW<I>, 

1-75,  160  >7 

88 

X 

1194, 

1157, 

1121, 

1086, 

1052, 

1020, 

89 

X 

9887, 

9583, 

9289, 

9004, 

8729, 

8463, 

88 

ft 

8206, 

7957, 

7716, 

7483, 

7257, 

7039, 

9:1 

X 

6828, 

6624 , 

6426, 

6235, 

6850, 

5870, 

92 

X 

5697, 

5529, 

5366, 

5209, 

5056, 

4908, 

93 

X 

4765, 

4626, 

4492, 

4262, 

4236, 

4114, 

94 

V,.' 

3995, 

3881, 

3769, 

3662, 

3557, 

3456, 

93 

X 

3358, 

3263, 

3170, 

3881, 

2994, 

2910, 

96 

X 

2829, 

2750, 

2673, 

2598, 

2526, 

2456, 

97 

X 

2388, 

2322, 

2258, 

2196, 

2136, 

2678, 

98 

X 

2021, 

1966, 

1912, 

I860, 

1810, 

1761, 

99 

X 

1714, 

1668, 

1623, 

1579, 

1537, 

1486, 

108 

X 

1456, 

1.418, 

1383, 

1349, 

131 6, 

1284, 

98 


1.01 

X 

1252, 

1.722, 

1192,  1162,  1134, 

1106, 

102 

X 

1079, 

10537 

103 

DATA  <[>N<  I  ), 

1-161,  229)7  1027  , 

9650  ,  96 

104 

X 

8056, 

7606, 

7179,  6774,  6390, 

6026, 

105 

X 

5681., 

5354, 

5044,  4751,  4474, 

4211, 

106 

X 

3963, 

3728, 

3505,  3295,  3097, 

2909, 

107 

X 

2732, 

2565, 

2407,  2258,  2118, 

1985, 

108 

X 

I860, 

1743, 

1632,  1528,  1430, 

1337, 

109 

X 

125G, 

1168, 

1092,  1019,  9517, 

8880, 

110 

X 

8283, 

7727, 

7197,  6704,  6237, 

5795, 

1:1.1. 

X 

5382, 

4997, 

4639,  4304,  3992, 

3702, 

11.2 

X 

3431, 

31.79, 

2945,  2727,  2524, 

2335, 

1 13 

X 

2160, 

1997, 

1.846,  1.705,  1575, 

1454, 

1:14 

X 

1342, 

1238, 

1141,  1052,  9694 

7 

115 

DATA  <DN<  I  ), 

1-  230,  302)7  8928, 

116 

X 

8220, 

7564, 

6958,  6366,  5824, 

5328, 

11.7 

X 

4875, 

4460, 

4081,  3734,  3416, 

3126, 

118 

X 

2860, 

261 6, 

2393,  2188,  2000, 

1828, 

119 

X 

1670, 

1526, 

1393,  1273,  1162, 

1061, 

120 

X 

9685, 

8842, 

8071,  7367,  6725, 

61 39, 

1.21. 

X 

5604, 

4695, 

3935,  3300,  2769, 

2325, 

127 

X 

1954, 

1643, 

1381,  1161,  9708, 

8111, 

173 

X 

6838, 

581 1, 

4975,  4289,  3720, 

3246, 

124 

X 

2847, 

2509, 

2222,  1977,  1767, 

1585, 

125 

X 

1428, 

1291, 

1171,  1065,  9717, 

8889, 

126 

X 

81  52, 

7494, 

6904,  6374,  5897, 

5465, 

127 

X 

5074, 

4719, 

4396,  4101,  3831, 

3584  7 

128 

DO  10 

KH=1, 5 

129 

IFCHI. 

GE.  H9TS<KH))G0  TO  20 

130  10  CONTINUE 

131  KH-KH-1 

1.32  20  Dl.H-HSTP  <  KH  > 

1.33  IDLH-DLH 

134  IHI;  <  I  NT  <HI7Dl  H)  )*IDLH 

1  :SS  J=  <  H I -HSTS  <  KH  )  > 7DLH+J0FS  <  KH  > 

136  30  IF< J.  LE.  276)00  TO  40 

137  TP2-TPC  J+l)71 0. 

138  IF< J.  EQ.  277) GO  TO  45 

139  TP1-  TP<  J)710. 

140  GO  TO  50 

141  40  TP2=TP<  J+D7100. 

142  45  TP1-TP  <  J  )  71.00. 

143  50  TNPW--5 

144  IF< J.  GT.  80 ) TNPW~-6 

145  1F(J.  GT.  161>TNPW~-7 

146  IF<  J.  GT.  198>TNPW-:— 8 

147  IF<  J.  GT.  228)TNPW- -9 

148  IF<  J.  GT.  254 ) TNPW--10 

149  IF<  J.  GT.  270)TNPW:--11 

150  IF<  J.  GT.  288>TNPW-  -12 


99 


1.51 

152 

ir.4 

1.55 
1 56 
137 
It  .8 


DN1=DN  <  J  >  *  <  10.  **TNPW> 

IF<  <  J.  EQ.  80>.  OR.  <J.  EQ.  161  >.  OR.  <J.  EQ.  198>>G0  TO  68 
IF<  <  J.  EQ.  228  >.  OR.  <J.  EQ..  254  >.  OR.  <J.  EQ.  2?0>>GO  TO  60 
IFCJ.  FGf.  288 >GO  TO  6® 

00  TO  70 

60  TNPU-TNPW— 1 

78  DN2---DN  <  J+l  >  *  ( 10.  **TNPW> 

RRT 1 G~  <  H I  —  IH I >  /DLH 


j.  3S# 
.1.60 
16:1 
J  62 
1 67 


T-  TP1 +RRT 1 0*  <  TP2-TP1 > 
RHG0=DN1* <  < DN2/DN1 >  **RRT 1 0 > 
100  CONTINUE 
RETURN 
END 


100 


1 

2 

♦ASSM 

SCRAT 

3 

BL.KDTA 

PROG  INITIALIZE  COEFFS 

14  NOV  79  R01 

S  REM3 :  BL.KDTA.  1 

4 

CROSS 

*> 

NORX3 

6 

Nl.STC 

7 

Nl.IST 

8 

$FORT 

9 

BLOCK  DATA 

.1.0 

COMMON  /C0FF/PXV1<51>,  PXV2<35>,  P22C21 >,  PZ1<  51 > 

11 

DATA  PXV1/ 

.  0049275309, 

±2 

X 

.  0032924339, 

.  0018334202, 

.  0005431455, 

-.  8005857183, 

13 

X 

-  0615605104, 

0023885688, 

-.  0030772295, 

©036338298, 

14 

X 

— .  0040657036, 

0043801926, 

0045846291, 

*.  0046863556, 

15 

X 

0046926998, 

064611.0079, 

©044486076, 

-.  0042128451, 

IS 

X 

0039110482, 

-.  0035565621, 

-.  0831387159, 

-.  8026828498, 

17 

X 

00219G3003, 

0016684636, 

-.  0011244968, 

-.  8005659165, 

18 

X 

.  0000000000, 

.  0605659165, 

.  0011244968, 

.  0016684036, 

19 

X 

.  0021903003, 

. 0626828498, 

.  0831387159, 

.  0835505621, 

28 

X 

.  0039110482, 

.  ©042128451, 

.  0044486876, 

.  0046110879, 

21 

X 

.  0046926998, 

. 0046863556, 

.  8045846291, 

.  8043801926, 

22 

X 

.  0040657036, 

.  0636338298, 

. 08307/2295, 

.  0023885688, 

25 

X 

.  0015605104, 

. 0005857183, 

-.  0085431 455, 

-.  0818334202, 

24 

X 

0032924339, 

-.  0049275309/ 

20 

DATA  PXV2/ 

.  009731 5572, 

.  0048963241, 

.  008836152 6, 

26 

X 

0024974160, 

0051528104, 

-.  0071784742, 

-.  0036228512, 

2? 

X 

0095343813, 

-.  0099615045, 

-.  0899526607, 

-.  0095562898, 

28 

X 

0088208355, 

-.  0077947341 , 

-.  0065264329, 

-.  0850643682, 

29 

X 

0034569885, 

-.  ©017527097, 

-.  0000000000, 

.  8017527097, 

30 

X 

.  0034569805, 

.  0050643682, 

.  0065264329, 

.  0077947341, 

31 

X 

.  0088208355, 

.  0095562898, 

.  0099526607, 

.  089961 5045, 

32 

X 

.  0095343813, 

.  0086228512, 

.  0871784742, 

.  0051528184, 

33 

X 

.  0024974160, 

-.  0008361526, 

-.  0848963241, 

-.  0897315572/ 

34 

DATA  PZ2/ 

.  0231351145, 

35 

X 

.  0027605407, 

0119105503, 

-.  0215118565, 

-.  0266771168, 

36 

X 

0280400403, 

-.  0262343623, 

-.  0218937919, 

01  56520642, 

37 

X 

0081428923, 

.  0000000000, 

.  0881428923, 

.  0156520642, 

38 

X 

.  0218937919, 

.  0262343623, 

.  0280400403, 

.  8266771168, 

39 

X 

.  0215118565, 

.  0119105503, 

-.  0027605407, 

-.  0231351145/ 

40 

END 

101 


$ASSM 


CONTRL 


SCRAT 

PROG 

CROSS 

N0RX2 

NLSTC 

NLISf 


CONTROLS  FLOW  30  NOV  79  R01  S  REM3 :  CONTRL.  FOR 


TFORT 


SUBROUTINE  CONTRL <GS, N2MID,  KSM2, KSN,  IPGF,  IV> 

COMMON  /CONST/  RG,  PHI,  HliSL,  BITS,  ZB,  DIR,  VB,  RM,  GSRG,  RB 
COMMON  /COE r'/  PXVK51  >, PXV2<25>,  PZ2<21>,  PZ1<51> 

COMMON  /THFO/  ZT<2>,  RZ<2>,  DZ,  XT<3>,  RX<2>,  DX,  VT<2>, 

X  RV<3>,  DV,  WX,  WV 

COMMON  /COR IS/  OMFGfi,  CPH,  CPHCZB,  CPHSZB,  CZB,  SPH,  SZ.B 
COMMON  T I ME < 1 00 > ,  X<100>,  V<100>,  Z<100>,  TMID<5©>,  XMID<50>,  VMIDC50 
R  ZM I D  i  58  ,  XVM  <  06  > ,  WM<5'3>,  ZVH<58>,  IX,  J,  N,  MIDMRX,  KMRX1,  KMRX2, 

B  CM,  NIST,  NKVi,  NZ1,  ZXM<5©> 

IHTEGFk+2  IX,  -J,  N,  MIDMRX,  KMRX1,  KMRX2,  NIST,  NXVl,  NZ1 
INTF.GER+2  I,  MM,  IV,  IPGF,  KSN,  KSW2,  N2MID 


D INFUSION  I GRR < 8©6 > ,  I COMSV < 866 > 

EQU I VALENCE  < T I ME < 1 > ,  I GRR  <1>) 

MM--NZ:l/2+i 

XT  <1>-X<'MI,1> 

VTCd  V  <  MM 
?T  < 1 > “Z < MM > 

XT  (  2  >  —  0.  0 
VT < 2 > “0.  0 
21  <  7.  >  -  0.  0 
00  1  l':  1,  NZi 

XT 2 >  =-XT < 2 > +PZ1  <  I  >*X<  I  >*?..  © 

VT  <  *■  >  r  VT  < 2  >  +PZ1  <  I  >  * V  <  I  >  *2.  0 
1  ZT < 2 > -ZT < 2 >  +PZd  <  I)*Z<I) *2.  0 
Q  S  R  Ci  - = G  S  /  R  G 
RM-  2.  0+AB/RMS 
Rflf>«57.  29527793 
CZB-  COS  •:  ZB/RRD  > 

SZB  ■  S I N  <  ZB/RRD  > 

OMEGA-  0.  0001458 
CPH- COS  < PH I /RAD > 

SPH  S  J  N  <  PH  I /Rfi£>> 

CPH  S  Z  B  ~  C  P  H  +■  S  2  B 
CPHCZB-  CPH*CZB 
WX-  6.  0 
WV---  0.  0 

DO  10  1-1,806 
10  I COMSV  <  I  >  -  I  GRR  <  I  > 

CALL  TROBINCGS,  N2MID,  KSW2,  KSN,  IPGF,  IV> 
DO  20  1-1,806 

20  i  grr  <:  i  :>  =■- 1  comsv  c  i  :> 

Rb TURN 


♦HSSM 


2 

SCRfiT 

3 

CORRS 

PROG  NUI1  OF  SMOOTH  PIS  14 

4 

CROSS 

5 

NGRX3 

6 

Nl.STC 

7 

NLIST 

8 

♦FORT 

9 

SUBROUTINE  CORRS <N,  ft,  11,  B,  C2,  IS> 

18 

DIMFNSION  8<1>,  B<1> 

±1 

INTEGER*2  I,  J,  K,  IN,  IR,  IW,  11,  N 

12 

D=0.  0 

13 

DO  1  I"l, N 

3  4 

1 

D  D+R<  I  >  =+-R  <  I  > 

IS 

C2--0.  © 

36 

DO  50  J- 1, 11 

17 

DO  50  K—l, n 

18 

IF < J.  ED.  K>  GO  TO  50 

19 

I R~ I BBS  <  I S*  <  J-K  > > 

20 

IN~N— IR 

23 

RI1SR"  0.  0 

22 

IFCIN.  LE.  0)  GOTO  5© 

23 

DO  10  I~l,  IN 

24 

IW---I  +  IR 

25 

10 

RMSR- RIISR+fK  I  >*ft<  IW> 

26 

RMSR“RI1SR/D 

27 

C2-C2+B  C  J  >  *B  C  K  >  *RI'1SR 

28 

50 

CONTINUE 

29 

RETURN 

.30 

♦flSSM 

31 

L  1ST 

32 

♦FORT 

33 

END 
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$ASSM 


1 


4 

5 

6 
7 

9 

aC-i 

ij 
1  2 
_i.  s' 

14 

ir. 

d  7 
18 
1.3 
2fi 
21 


24 

23 


29 
7.0 
70 
72 
S  s 
74 
33 

36 

37 

38 
7;9 
48 

41 

42 

43 

44 
43 


SCRRT 

DFCAl.T  PROG  SPHFRF  DROPPING  14  NOV  79  R61  5  RfcM3 :  DECALT.  FOR 

CROSS 
N0RX3 
NLSTC 
Nl  1ST 

$FOR  l 

SUBRGUTINF  DFCAl.T ( JCNT,  IFLAG1> 

COMMON  T I MF (.  1 60 > X <  1  00 > ,  Y<100>,  ?<100>,  TMir><50>,  XMID<50>,  VMID<50>, 
0  ZMItX30>,  XVM<5©>>  VVM<50>,  ZVM<50>,  IX,  J,  H,  MIDMAX,  KMRX1,  KMRX2, 

B  Cl,  NIST,  NX VI,  NZ1,  ZXM<58> 

I NTEGF 8*2  IBKD,  ICNT,  JCNT,  IFLAG1 

INTEGER*2  IX,  J,  N,  MIDMRX,  KMtiXl,  KMRX2,  NIST,  NXV1,  NZ1 
I FLAG!  -1 

R-  6271229.  313 
SI  1-  8. 

S71-C 

S31--0. 

C  DATA  READ  Six-X,  S2~V,  S3«2 

1  RFRDC1 ,  20,  FND-"8>  T,  SI,  S2,  S3 
20  FORMAT C4F18.  7> 

D I F  S-  <  <  S  J  -SI  1  >  *  <  Sl+Sl  1  >  +  <  S2-S21  >  *  <  S2+S21 !)  > 

DOTH  S31-S2-0 1 FSX <  R*6.  361 666  > 

IFCDLTH>  2,2,3 

2  ICNT  -6 
G  TO  4 

2  ICN i  ICNT+1 

IF <3- ICNT >3,  3,  4 
4  Six  SI 

S3!  ;  S3 
GO  TO  1 

3  I  AND-- 7 
TM-T+10.  0 

FT--  <  IFIXC TM >  + 1 BND > /10. 

6  READ  (!  1,  20,  END- 8 >  T,  SI,  S2,  S3 
IFCFT  -  T>  7,7,6 

7  RETURN 

8  IFLAG1-0 
RETURN 

1ASSM 

LIST 

SFORT 

END 


1  04 


*HSSM 


rf 


y 


1 

2 


< 

4 

5 


8 

9 

10 
11 
1,? 
:1  3 

14 


DFV 


SCRfiT 

PROG  POS,  DENSE,  &WIND  CHG 

CROSS 

N0RX3 

NL  STC 

NL  1ST 


14  NOV  79  R01 


RE  M3  :  DEV.  FOR 


$FORT 


SUBROUTINE  DEVCVRHG,  VWX, VWV,  VP,  RHOOG> 

COMMON  T I ME < 1 00 > , X  <1.005, V<108>, Z<100>, TMID<50>,  XMID<50>,  VMIDC50 
ft  ZM I D  <  50  > ,  XVM  <  50  > ,  VVM  <  50  > ,  ZVM<50>,  IX,  J,  N,  MIDMftX,  KMflXl,  KMAX2, 

B  Cl,  NIST,  NXV1,  N21,  ZXMC58) 

COMMON  ZD/X1,  X2,  X3,  VI,  V2,  V3,  ?1,  Z2,  Z2X,  Z2,  Z2X,  CONX1,  CO  NX  2,  C0NX3, 
ft  C0NX4,  CONZ1,  C0NZ2,  CGN73,  C0NZ4,  GS,  IDEGX2,  IDEGX1,  IDEGZ2,  IDF.GZ1 


1  5 

E 

N2MID 

1 6 

INTFGER+2  IX,  J,  N,  MIDMftX,  KMflXl,  KMftX: 

1  7 

INTFGFR+-2  IV,  IDEGX1,  IDEGX2,  IDEGZ1, 

18 

DftTft  IV/0/ 

1  9 

VftRR;-36 

20 

VflRE  :  15 > * < .  15 >  +.  96383E-6 

21 

VftRft-VftRE 

22 

TWOSQ-  XI *X1+V1*V1 

23 

THRSQ--  TWOSQ+Z1*?! 

24 

VftRX-  X1*X1*VFIRR/THRSQ+X1*X1*21*Z:]  .+' 

■-.it  ( 

VflRV-:  VI  *Vl*VftRR/THRSQ+Vl  *V1  *Z1*Z1*' 

26 

V  ft  R  ZZ 1  *  2 1  *  V  ft  R  R T  H  R  S  Q + T  W  0  S  Q  *•  V  ft  R  F 

i** 

I F  <  I DEGX1.  NF.  1>  GO  TO  31© 

28 

VflRXV- CONXl+VflRX 

29 

V  ft  R  V  V  -C  0  N  X 1 + V  ft  R  V 

30 

V  ft  R  X  7 V = C 0 N XI* V ft R Z 

31 

GO  TO  32© 

.  S  C.' 

31  0 

V  ft  R  X  V  ;C  O  N  X  2  *  V  ft  R  X 

”s": 

V  ft  R  V  V  -C  0  N  X  2  *■  V  ft  R  V 

.34 

VftRXZV*  C0NX2+VftRZ 

35 

320 

I F  a  DEGX2.  NF.  1>  GO  TO  33© 

36 

V  ft  R  X  ft  -  C  0  N  X  3  ■+:  V  ft  R  X  V 

r 

V  ft  R  V  ft-C  0  N  X  3  *  V  ft  R  V  V 

V  ft  R  X  Z  ft  -C  0  N  X  3  *  V  ft  R  X  2  V 

“'.9 

GO  TO  340 

401 

330 

V  ft  R  X  ft--  C  0  N  X  4  *  V  ft  R  X  V 

41 

V  ft  R  V  ft  ;C  0  N  X  4 + V  ft  R  V  V 

43 

VfiRX2ft~CGNX4*VfiRXZV 

43 

340 

IFCIDFGZl.  NF.  1>  GO  TO  350 

44 

VftRZV-SQRT < COUZl+VftRZ > 

45 

GO  TO  360 

46 

35© 

VftRZV-SQRT < C0NZ2+VARZ > 

47 

360 

IF<  IDEGZ2.  NF.  1>  GO  TO  370 

48 

VflRZft- SORT  <  C0N23*VflRZV**2  > 

49 

GO  TO  380 

50 

37© 

VftRZft— SORT  C0N24*VfiRZV**2  > 

N2MID 


105 


4.>£C*iGSP 


51  sm  VRHO— <  2.  *VRRZV/Z2  >  **2+  <  VRRZRX  <  Z3-GS  >  > **2 

52  IF  CIV.  EQ.  1>  GO  TO  390 

03  VP— GS**2X4.  6*VRH0*  <  ZM I D  <  N2M I D  >  -ZM I D  <  N2M I  D+l  >  >  **2*RH000**2 


54  IV-1 

55  GO  TO  40O 

h€  390  VP- VP+GS++2/4 .  0* VRHO*  C ZM I D  <  NPM I  D-i  > -ZM I D  < N2M I D+l >  > **2*RH000**2 

5?  400  VRHO:- SORT C VRHO > 

.1?  VWX=  SORT  <  VRRXV+  <  Z2X/  < Z3X-GS  >  > **2* VRRXR+  < X2X <  Z3X-GS  >  )  **2*VRRXZV+ 

59  1 C  X3+-Z2X/ < Z3X-GS > **2 > ++P+VRRXZR  > 

m  VWV-  SORT  <  VRRVV+  <  Z2X/ <  23X-GS  >  > **2+Vfl RVR+  < V3/ <  Z3X-GS  >  > **2*VRRXZV+ 

61  1  <  V3+Z2X/  <  Z3X-GS  >  **2  '>  **2*VRRXZfl  > 

62  RETURN 

63  $RSSM 

64  I.  1ST 

65  If-  OPT 

66  END 
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1 

2 

3 

4 

5 
€ 
f 

8 

9 

10 

1:1 

:1.2 

13 

14 

15 

16 
1? 
18 

19 

20 
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$RSSM 

SCRAT 

DRRGT  PROG  DRAG  TABLF  14  NOV  79  R01  S  REM3 :  DRAGT.  FOR 
CROSS 
N0RX3 
Nt.STC 
Nl  1ST 
*FORT 

SUBROUTINE  DRRGT CCD,  AMC,  RF,  I FLAG > 

C  ROUT  INF  SENT  RE,  AMC  :  RF.VNOL.DS  NO.,  MACH  NO. 

C  RETURNS  CD  :  DRAG  COEF. 

C  I FLAG-  -1  VALID  CD;  I FLAG  -  O  NO  CD 
DIMENSION  DRAG <2,  2> 

INTE.GER+2  AR1C24,  22  >,  AR2C17,  31  >,  AR3C16,  15  > 

INTEGERS  RFNL. C 40 > ,  AMACHC53> 

INTFGFR*2  ARIA C 24,  8>,  AR1BC24,  7>,  AR1CC24,  7) 

INTEGERS  AR2AC17,  S),  AR2BC17,  7>,  AR2CC17,  8>,  AR2DC17,  8) 
INTEGERS  I,  J,  16,  J14,  II,  JJ,  NMR,  I  FI  AG,  II,  J1 
EQUIVALENCE  CAR1C1,  1>,  AR1AC1,  1>>,  CAR1C1,  9), AR1BC1,  1 >>, 

X  CAR! Cl,  1.6  >,  AR1CC1,  1>  >,  CAR2C1,  1>,  AR2AC1,  1>  >, 

X  CAR2C1,  9>,  AR8BC1,  1>>,  C  AR2C1,  16>,  AR2CC1,  1>  >, 

X  CAR2C1,  24>,  AR2DC1,  1>> 

C 

DATA  ARIA/ 

X1673,  1604,  1535,  1491,  1449,  1384,  1342,  1289,  1250,  1199,  1171, 
X1132,  1104,  1075,  1069,  1055,  1038,  0,  0,  0,  0,  0, 

X1633,  1612,  1541, 1491,  1452,  1392,  1347,  1292,  1252,  1200,  1171, 
X1134,  1.106, 1078,  1069,  1055,  1041,  0,  0,  0,  0,  0, 

X1697,  1622,  1547,  1495,  1458,  1400,  1352,  1298,  1254,  1202,  1171, 
X1136,  1108,  1080,  1067,  1055,  1043,  0,  0,  0,  0,  0, 

X1720,  1640,  1555,  1503,  1463,  1410,  1358,  1302,  1259,  1207,  1175, 
X11.38,  1110,  1082,  1065,  1055,  1045,  0,  0,  0,  0,  0, 

X1749,  1658,  1567,  1519,  1473,  1420,  1366,  1316,  1266,  1212,  1184, 
X1141,  1111,  1088,  1070,  1058,  1051,  0,  0,  0,  0,  0, 

X1780,  1680,  1580,  1533,  1485,  1430,  1375,  1319,  1275,  1220,  1192, 
X1149,  1122,  1095,  1076,  1062,  1055,  0,  0,  0,  0,  0, 

X1S09,  1703,  1597,  1548,  1498,  1444,  1390,  1330,  1287,  1230,  1202, 
X1158,  11.30, 1101, 1084, 1066,  I860,  0,  0,  0,  0,  0, 

XI 841,  1728,  1615,  1564,  1.512,  1456,  1400,  1341,  1298,  1242,  1214, 
X1170,  1140,  1109,  1097,  1084,  1G72,  0,  0,  0,  0,  0, 

DATA  AR:1F:/ 

X1870,  1753,  1636,  1582,  1528,  1471,  1414,  1355,  1310,  1257,  1227, 
X1179,  1.146,  1113,  1100,  1089,  1075,  1063,  1050,  0,  0,  0, 

XI  903,  1780,  1657,  1601,  1545,  1487,  1429,  1368,  1322,  1270,  1237, 

X1186,  1.1.53, 1120,  1106, 1092,  1079,  1065,  1652,  0,  0,  0, 

X1941,  1810,  1679, 1621,  1562,  1502, 1442,  1382,  1.337,  1282,  12^6, 
X1192, 1157,  11.22,  1108,  1094,  1080,  1066,  1053,  0,  0,  0, 

X1.958,  1825,  1692,  1632,  1572,  1511,  1450,  1390, 1342,  1288,  1251, 
X1195,  1160,  1124,  1110,  1096,  1081,  1065,  1053,  1046,  1038,  1025, 
X1976,  1841,  1706,  1644,  1582,  1520,  1458,  1397,  1350,  1293,  1255, 


1142, 

0, 

1144, 

0, 

1147, 

0, 

1150, 

0, 

1157, 

0, 

1163, 

0, 

1173, 

0, 

1186, 

0/ 

1196, 

0, 

1203, 

0, 

121.0, 

0, 

1213, 

1015, 

1217, 
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XI 199,  .1:162,  1125,  1111, 
X1999,  1860,  1721,  1657, 
X:1 200,  11 63,  Id  26,  11 11, 
X2023,  1880,  1737,  1671, 
X1202,  1163,  1125,  1110, 

DRTR  RR.1C/ 

X2046,  1901,  1736,  1687, 
X1203,  1164,  1124,  11.09, 
X2077,  1927,  1777,  1705, 
X 1.202,  1162,  11  22,  1108, 
X209O,  1940,  1796,  1717, 
XI 2  0G,  1158,  1120,  1105, 
X2107,  1956,  1805,  1731, 
XI 199,  1.157,  111.9,  1103, 
X2123,  1972,  1821.,  1746, 
XI 196,  1.154,  1116,  1102, 
X2142,  1990,  1338,  1762, 
XI 134,  1.152,  1112,  1096, 
X2169,  20G7,  1859,  1780, 
X 1.1 91,  1151,  1110,  1095, 

DRTR  RR2TX 

X1442,  1380,  1297,  1245, 
XI 066,  1031,  1036,  1027, 
X1443,  1378,  1292,  1240, 
XI 061 , 1046,  1031,  1G22, 
XI 4  43,  1  j-76,  128/,  1 2  4, 
X1054,  1041,  1034,  1018, 
XI 4  44 ,  1373,  1280,  1227, 
X1619  1032,  1026,  1010, 
XI- 45,  1369,  1273,  1218, 
Xl.0-14,  1027,  1016,  1803, 
X1445,  1365,  1266,  1269, 
XI 637,  1020,  1011,  996, 
X1444,  1339,  1257,  1193, 
X1027,  101 1,  1001 ,  989, 
XI 4 42,  1355,  1249,  1186, 
X101.6,  1000,  991,  977, 

DRTR  RR2B/ 

XI 440,  1347,  1238,  1172, 
XI 062,  936,  980,  968, 
X1437,  1339,  1230,  1160, 
X  987,  971,  960,  948, 
X1436,  1.327,  1214,  11.43, 
X  964,  954,  942,  929, 
X1432,  1315,  1197, 1125, 
X  947,  931,  919,  906, 
X1426,  1298,  1130,  1106, 
X  936,  918,  907,  892, 
X1414,  1231,  1156,  1080, 


1097, 

1082, 

1064, 

1054, 

1592, 

1529, 

1465, 

1483, 

1096, 

1082, 

1064, 

1052, 

1604, 

1533, 

1472, 

1410, 

1095, 

1081, 

1062, 

1651, 

1618, 

1549, 

I486, 

1418, 

1094, 

1080, 

1062, 

1050, 

1633, 

1562, 

1490, 

1422, 

1092, 

1078, 

1063, 

1048, 

1644, 

1568, 

1496, 

1427, 

1087, 

1074, 

1061, 

1046, 

1658, 

1580, 

1562, 

1420, 

1888, 

1073, 

1058, 

1042, 

1671, 

1590, 

1510, 

1432, 

1886, 

1070, 

1054, 

1039, 

1686, 

1603, 

1519, 

1425, 

1081, 

1066, 

1051, 

1026, 

1  700, 

1613, 

1527, 

1.438, 

1079, 

1062,  1047, 

1031, 

1210, 

1179, 

1141, 

1106, 

1020, 

1010, 

997, 

987, 

1205, 

1172, 

1127, 

1100, 

1617, 

1005, 

991, 

981, 

1197, 

1167, 

1122, 

1095, 

1011, 

1000, 

985, 

975, 

1190, 

1160, 

1120, 

1090, 

1004, 

9 92, 

978, 

967, 

1184, 

1154, 

1117, 

1086, 

992, 

983, 

971, 

968, 

1174, 

1146, 

1168, 

1079, 

985, 

977, 

962, 

951, 

11.62, 

1136, 

1898, 

1072, 

977, 

969, 

952, 

942, 

1150, 

1123, 

1081, 

1061, 

967, 

958, 

940, 

929/ 

1132, 

1187, 

1868, 

1046, 

956, 

945, 

928, 

917, 

1115, 

1085, 

1050,  1030, 

940, 

929, 

918, 

897, 

1G95, 

1070, 

1024, 

1011, 

924, 

910, 

892, 

878, 

1076, 

1052, 

1014, 

998, 

906, 

888, 

867, 

854, 

1057, 

1032, 

993, 

968, 

881, 

870, 

855, 

841, 

1027, 

1013, 

972, 

946, 

1047,  1039,  1026,  1016, 
1356,  1298,  1259,  1219, 
1046,  1027,  1025,  1016, 
1361,  1302,  1262,  1221, 
1044,  1034,  1022,  10122 

1368,  1306, 1265,  1223, 
1042,  1032,  1019, 1009, 
1372,  1307,  1265, 1222, 
1042,  1022,  1020,  1010, 
1374,  1306,  1264,  1221, 
1029,  1023,  1017,  1007, 
1376,  1205,  1262,  1220, 
1036,  1023,  1014,  1004, 

1278,  1302,  1260,  1217, 

1022,  1022,  1010,  1000, 

1279,  1200,  1258,  1215, 
1028,  1019,  1006,  996, 
1379,  1298,  1252,  1211, 
1023,1014,1001,  991 X 


1082, 

1078, 

1072, 

1865, 

1061, 

1053, 

1044, 

1024, 


1021, 

1005, 

986, 

968, 

947, 

926, 


301 

X 

910, 

895, 

882, 

864, 

854, 

846, 

828, 

813, 

102 

X1277, 

1247, 

13.19,  1052, 

1008, 

987, 

945, 

936, 

891  j 

3  02 

X 

869, 

857, 

843., 

829, 

818, 

868, 

796, 

783/ 

1B4 

DRTR 

RR2C/ 

3.05 

X1175, 

1062, 

952, 

900, 

857, 

837, 

798, 

773, 

750. 

106 

X 

727, 

715, 

705, 

696, 

686, 

677, 

667, 

656, 

107 

X1088, 

984, 

880, 

825, 

787, 

770, 

721, 

781, 

678. 

108 

X 

661, 

652, 

628, 

624, 

620, 

610, 

662, 

593, 

109 

X1027, 

928, 

825, 

770, 

727, 

710, 

677, 

650, 

625 

110 

X 

613., 

600, 

589, 

579, 

575, 

564, 

558, 

551, 

311 

X 

996, 

899, 

784, 

731, 

692, 

670, 

625, 

607, 

584. 

13.2 

X 

571, 

562, 

553, 

541, 

535, 

530, 

525, 

519, 

132 

X 

974, 

874, 

752, 

700, 

657, 

627, 

682, 

571, 

549, 

114 

X 

53?, 

526, 

53.8, 

513, 

505, 

500, 

499, 

494, 

13.5 

X 

955, 

852, 

728, 

674, 

622, 

614, 

575, 

542, 

527. 

136 

X 

535, 

504, 

498, 

488, 

482, 

479, 

473, 

470, 

117 

X 

938, 

824, 

708, 

.657, 

63  3, 

588, 

555, 

529, 

515. 

118 

X 

498, 

486, 

481, 

47?, 

474, 

469, 

462, 

457, 

119 

X 

925, 

819, 

689, 

639, 

598, 

572, 

54Q, 

517, 

500. 

128 

X 

489, 

479, 

470, 

463, 

460, 

456, 

450, 

448/ 

12:1 

DRTR 

RR2CV 

122 

X 

914, 

810, 

683, 

630, 

588, 

561, 

527, 

507, 

492. 

122 

X 

480, 

472, 

459, 

452, 

456, 

445, 

442, 

438, 

124 

X 

905, 

798, 

672, 

620, 

580, 

553, 

516, 

498, 

484. 

125 

X 

472, 

464, 

452, 

443, 

442, 

436, 

434, 

430, 

126 

X 

884, 

789, 

664, 

606, 

572, 

546, 

510, 

491, 

475. 

127 

X 

465, 

456, 

445, 

436, 

435, 

426, 

425, 

421, 

128 

X 

874, 

783, 

659, 

599, 

564, 

548, 

502, 

484, 

467. 

129 

X 

45?, 

448, 

427, 

429, 

426, 

43  9, 

436, 

411, 

128 

X 

85?, 

778, 

652, 

590, 

558, 

529, 

494, 

472, 

454 

123 

X 

445, 

425, 

428, 

421, 

434, 

407, 

402, 

395, 

122 

X 

849, 

762, 

649, 

584, 

543, 

518, 

486, 

461, 

444. 

122 

X 

425, 

424, 

419, 

43.0, 

405, 

399, 

393, 

384, 

124 

X 

846, 

758, 

648, 

577, 

525, 

505, 

470, 

449, 

432. 

125 

X 

425, 

416, 

406, 

298, 

393, 

291, 

384, 

377, 

126 

X 

844, 

751, 

647, 

575, 

530, 

492, 

460, 

444, 

430, 

127 

X 

414, 

407, 

396, 

389, 

387, 

384, 

380, 

372/ 

128 

C 

129 

DRTR 

RR3/ 

140 

X584,  575, 

0,  0, 

0, 

0, 

0, 

0,  0, 

0, 

0, 

141 

X544 ,  527, 

0,  0, 

0, 

0, 

0, 

0,  0, 

0, 

0, 

142 

X512,  507,  503 , 500, 

499, 

0, 

0, 

0,  0, 

0, 

0, 

142 

X490,  486,  483 , 480,  480, 

0, 

0, 

0,  0, 

0, 

0, 

144 

145 

146 

147 

148 

149 
158 


0, 

0, 

0, 

0, 


0, 

0, 

0, 

0. 


0, 

0, 

0, 

0, 


X464,  458,  457,  458,  458,  459,  461,  465,  469,  475,  477, 
X45i,  449,  446,  449,  450,  451,  454,  456,  459,  462,  468, 
X444,  442,  442,  441,  442,  442,  444,  447,  450,  451,  462, 
X437,  422,  423 , 422,  422,  424,  428,  440,  442,  442,  456, 
X429,  424,  425,  426,  427,  428,  421,  422,  424,  425,  449, 
X420,  41 6,  420,  422,  425,  426,  428,  421,  422,  422,  444, 
X410,  408,  416,  420,  422,  424,  425,  426,  427,  421,  442, 


485,  497,  502, 
482,  487,  495, 
472,  479,  485, 
464,  472,  479, 
456,  465,  474, 
450,  459,  468, 
448,  455,  462, 


0, 

0, 

0, 

0, 

513, 

501, 

491, 

485, 

480, 

475, 

468, 


0, 

0, 

0, 

0, 

512, 

505, 

496, 

492, 

486, 

480, 

475, 
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1  51  X394.  39  5,  .402.  40?.  410.  41 2.  414.  417.  420.  479.  435.  441.  445.  454.  45?.  465. 

I*-'.'  X_s.'8_  -  762,  391.  396.  398.  400.  406.  408.  410.  41.?.  428.  436.  440.  449.  451.  456. 

153  X376,  376.  385.  391.  393.  395.  400.  407.  404.  40?.  41 9.  431 . 439.  448.  450.  455. 

154  X371  •  370.  379.  386.  389.  391.  394.  396.  399.  407.  410.  425.  438.  44?.  450.  435/ 

1 *.  .5  C 

156  ORTR  RFNI./ 

1‘*.7  >;  i  2,  3.  4.  5.  -1.  10.  15.  20.  30.  40.  50.  60.  80.100.120. 

1'  .8  vJ  *>  '  j.60.  180.  200-  270.  250.  308.  350.  400.  450.  500.  550.  600.  650.  ?00.  800. 

1  7  *  y.  80>3.  1000.  1500.  2000.  2500.  3000.  3500.  4000/ 

J  60  C 

J  61  DHTC<  rlHRO  ' 

J  r.2  :  1.  'i  140.  .170.  400.  380.  360.  340.  376.  300.  280.  260.  250.  240.  230.  220.  216. 

a -74  *10  1  1 90.  185.  180.  1  75,  176.  165.  166.  155.  150,  145,  140,  135,  130,  125, 

164  :  1  0  ‘  V ..  xl0,  165,  100,  95,  98,  80.  30.  75.  70,  65,  60,  55,  50,  45, 

165  X  40  ■  !•  70,  10,  0/ 

]  26  f 

167  C 

1  -  °  if  1  Fr  j"T  j. '  i0.  6  OR  RHC.  GT.  1  81  GOTO  10 

I-"’1*  ! 

i  ,'ij  GOTO  70 

1,1  16  I F  ■  RF  -16  :;0.  30,  36 

1 72  20  1=1 

1  i  3  GOT  14  70 

X 74  30  IF'  RF -4 0000  )  50,50,40 

175  40  I • 77 

176  GOT ' 1  76 

1  77  50  f  '1  60  16  1  ■  -  0 

178  I-  16 

179  TR1  PrNLC  I  +1  >*10. 

180  IF'  TP  i  i  0.  TR1*75. 

18:1  IF'.  Fir  Ir  i  >  70,60,60 

187  60  CONTINUE 

185  70  IF  . Hl'IC-4.  S'  90,90,80 

184  80  J=  1 

185.  GOTO  110 

186  90  00  j.OO  J141.57 

187  J--.T14 

188  Tkl  RHRCHC  J+D/100. 

189  IF  t.  Rl  IC-TR1  >  100,  1 10,  11  6 

190  1610  CONTINUE 

.1 91  C 

197  C  DRRG  TRBLE  SPLIT  INTO  THREE  SUR-RRRRVS 

193  C  Ri-:1,  RR2,  RR3 

194  C  IF  NOT  IN  OF  THESE •  COEF=-0 

195  C  FIND  FOUR  COEK:  Cl,  J>,  <  1+1,  J>,  <1,  J+l  >,  <  1 4-1,  J+l> 

196  C 

197  110  DO  200  ri:- 1, 2 

198  DO  200  N-  1 ,  2 

1.99  II=I+M-1 

7  OH  T.J--J+N-1 


no 


'1 


229 

230 

231 

232 

233 

234 
23'. 
23* 
23? 
238 
233 

240 

241 
24? 

243 

244 

245 
24* 

247 

248 


201 

NMfl-4 

202 

IF<  JI.  LE  24)  GO  TO  120 

203 

IF<  JJ.  LT  39)  GO  TO  140 

204 

NMR=  3 

20b 

GO  TO  14© 

20* 

120 

IF< JJ.  GT  22)  GO  TO  130 

207 

NMR=1 

208 

00  TO  140 

209 

120 

IF<  1 1.  LT.  8)  GO  TO  140 

210 

NMR  -2 

211 

140 

GO  TO  <150,  160,  170,  220) 

212 

150 

DGV  -'RRI  <11,  J  J  )  31000. 

23  3 

GO  TO  190 

2:1  4 

160 

13=  IT-7 

21  5 

Jl=-  JJ— 22 

21  * 

OGV  RR2  <  11,  Jl)  21000, 

217 

GO  TO  190 

2:18 

170 

11=  11-24 

219 

Jl~  JJ-38 

220 

DGV  RR3<  11,  J3  )  21000. 

22.1 

190 

CONTINUE 

t '*  3' 

DRRG < M,  N) - DGV 

22.3 

200 

CONTINUE 

2  24 

C 

225 

TR1=  RFNL<I)*10 

22* 

IF<TR3  LT  0  >  Tkl-75. 

227 

Tk2-RFNL<I+l)*10. 

228 

IF<TR2  LT  0.  )  TR2-75. 

NMR 


230 


22© 


500 
♦RSSM 
L  1ST 
♦FORT 


Rl=  <  TR1-RF  )  2  <  TR1-TR2  ) 

Tkl=  RMRCH<J)210© 

TR2=^RMRCH<  J+l 1 /100. 

R2  -  <  TR1-RMC >2 < TR1-TR2  > 

IF<DkRG<2,  1)  EQ  0.  0.  OR.  DRRG<3,  1>.  EQ.  0.  0  OR  DRRG<2,  2>.  EQ  0.  0 
X  .  OR  DRRG<1,  2)  EG!  0.  0)  GOTO  220 
CD1-R1*  <  DRRG  <  2,  1)-DRRG<1,  1>>+DRRG<1,  1) 

CD2-R1* < DRRG < 2,  2)-DRRG<l,  2>)+DRRG<l,  2) 

CD-R2*  <  CD2-CD1 > +CD1 

IFLRG  -1 

RETURN 

1RE-RF 

WRITF<4,  500)  RMC,  IRE 

IFLRG  0 

RETURN 

F0RMRT<238X,  20HDRRG  VRLUE  NOT  GJVFN,  14X,  F5.  2,16) 


END 
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♦HSSM 


DRVT 


^  FORT 


SCRAT 

PROG  FINDS  DRAG  VALUE  14  NOV  79  R03  S  RBM3.DRVT.  FOR 

CROSS 
N0RX3 
NLSTC 
Nl  1ST 

SUBROUTINE  DRVTCK,  TEMP,  RE,  AMC,  CD,  DENTT,  HI,  LLL) 

COMMON  /CONST/  RG,  PHI,  HMSl,  AMS,  2H,  DI,  VB,  AM,  GSRG,  AB 
COMMON  /COR IS/  OMEGA,  CPH,  CPHC2B,  CPHSZB,  CZB,  SPH,  S7.B 
COMMON  ,'THF.O/  20),  RZ<2>,  DZ,  XC3),  RX<3>,  DX,  VO),  RV<3>,  DV,  WX,  WV 
COMMON  /T r  r:l  /  I  MX,  B I  ASH  I  <  600  ) , 

•<  h  I  ASM/  <  600  ) ,  b  I ASMV  <  €06  ) ,  B I ASR  <  600  ) ,  B I  AST  <  600 ) ,  JMX,  B I ASMC  < 600 ) 
I  "J ;  i  GEP+2  K,  LI  L.  L,  I  Ft.  AG 
DATA  L  ■  OtUT/2.  0.  0/ 

H I  -  R?  <  V  >  +  <  RX  <  1 )  *RX  (  1  >  +RV  <  3  >  *RV  <l>>/<2.  *RG  )  ♦HMSL. 

OP- 1.  +HI/RG 

WT HR  GSRG/ ( nP**3 > 

WTWO- t J  .  -VB+ DENT/AMS > ♦WTHR 

BX'-'PX  '•  j  .♦Wind 

BV"  R V  J  />*MTMO 

B2  R21-.  1  '  ►  G  ) +T  1 T  WO 

CX=  OMEGA + • RT < 7 > *  CPHS2B+RV < 7 > ♦SPH  > 

C V '  OMF OR  R2 1 1  ’>  +  CPHCZB— RX  <.  7. )  ♦SPH  ) 

02= -OMF  GA*CPH+ < RX <2> ♦SZB+RV <2> ♦CZB ) 

IF*  *  FO  2)  GOTO  2 

VEL  SQR T < RZ ■  2 > ♦RZ C 2 >  +  < RX < 2 ) -  WX > ♦♦2+  <RV<2> -WV > ♦>*<2 > 

C  HI  L  AT MOS (Ml ,  TF MP,  DENT  ) 

RE-  f'F !'•!'.  +  VFI .  +-D I  *  <  Tfc MP+110.  4>/<3.  458F-6*TEMP**i.  5> 

AMC  VFL  +  SORT  (288.  1 6/TEMP )/340.  29205 
CALL.  DPAGT  < CD,  AMC,  RE,  I  FLAG 
IF1  I.  1  AO  200 ■  1G6,  206 


10«?| 

NR i  Th  < 4  'TO  :■  HI 

25 

90© 

FORMAi <1X,  19HWARNING 

26 

/ 

i  3X, 37HZER0  IN 

37 

2wGi 

com  ini. if 

/R 

E  0  5  ♦  D E.  N  T + C  D ♦  AM  ♦  VFL 

39 

DZ  RZ'  2> 

40 

D  f  • :  ( 2  :> 

43. 

DV-  RV<?> 

42 

GOTO  6 

43 

DZ:--E*R?  <  7  )  +BZ-CZ 

•14 

DX:  -F+ < RX ( 2 ) -WX > +BX- 

45 

DV  -  -E*  (  R V  (  7  >  — W’ r‘ )  +BV- 

46 

2</> ~ DZ 

47 

X < 2 > :  DX 

48 

V  <  .2  >  ••  DV 

49 

LLL  l 

‘  10 

DENT  M- DENT 

51 

6 

RETURN 

r  * 

*«SSI1 

r  •  “• 

LIST 

54 

*F0R1 

r  •  r~ 

.  i  \ 

END 
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♦ASSN 


1 


4 

5 

6 

7 

8 
9 

10 

11. 

18 

13 

14 
iu 
16 
17 
1  8 
19 
80 
81 
88 
23 
84 
8l*> 
86 
8? 
88 
89 

30 

31 
38 

33 

34 
30 

36 

37 

38 

39 

40 

41 
48 

43 

44 

45 

46 

47 

48 

49 

50 

51 
08 
03 
04 


SCRAT 

FI  TON  PROG  CALL  SMOOTHFD  VEl  14  NOV  79  R01  S  REM3  :  F I  TON.  FOR 

CROSS 
N0RX3 
Nl.S'l  C 
NLIST 

♦FORI 

SUBROUTINE  FI TON 

COMMON  /COEf- /  PXV1  < 0 1  > ,  FXV2 <3G>,  P22 <  21  > ,  P21  < 30 > 

COMMON  TIME <100 1.-  X<100>,  V<10G>,  Z<100>,  TMID<S0>,  XMID<50>,  VMID<50>, 
ft  2M I D < 50 > ,  XVM<50>,  VVM<50>,  ZVM<50>,  IX,  J,  N,  MI  DMAX,  KMAX1,  KMAX2, 

B  Cl ,  N1 ST,  NXV1,  NZ1,  7XM<50> 

INTEGER*?  IX,  J,  N,  MIDMAX,  KMAX1,  KMAX2,  N1S  f,  NX VI,  N71 

INTEGER*?  I,  K,  NFND 

PX=;  0.  0 

PV=0.  0 

PZX=  O. 

P&  0.  O 

DO  1  I  I,  KMflXl 
PZX^ P2X+2< I >*PXV1< I ) 

PX~PX+X< I >*PXV1< I > 

1  PV  PV+V  < I >  *PXV1 < I > 

NFND-  NZ1+NJ  ST-1 

DO  2  I  •••NIST,  NEND 
K>  I-N1ST+1 

2  PZ-PZ+2< I >*PZ1 <K> 

GO  TO  5 

ENTRV  FI TON? 

PX-0.  0 
PZX’  0. 

PV-  0.  G 
PZ:  0.  0 

DO  3  L  I,  KMAX1 

3  PZ=-P2+2<  I  >*PZ1  <  I  > 

NEND-NXV1+N1ST-1 
DO  4  I -NIST,  NFND 
K::  I-N1ST+1 

PZK  PZX+Z  <  I  >  *PXV1  <  K.  > 

PX=  PX+X  < I >  +PXV1 <K> 

4  PV- PV+V < I > *PX VI <K> 

5  7MID<N1=TIME <MIDMAX) 

XM I D  <  N  >  —X  <  M I DMAX 1 
VMID<N>~V<MIDMRX> 

ZMID<N>~Z<MI DMAX > 

XVM<N>~PX*8. 

VVM  <  N  >  :rPV*2. 

2 VM  <  N  >  -  PZ*2. 

ZXM<N>~PZX*2. 

RETURN 

♦ASSM 

LIST 

♦FORT 


END 
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#RSRM 


INTER 


i  P  UR  I 


scrrt 

PROG  DRTR  I NTERPOLRT I  ON  14  NOV  79  R01  S  RFM3 : INTER.  FOR 

CROSS 

nor::  s 

Nl.STC 

NLIST 

SUOROIJT I NP  I N  l  ER  <  I NTERR,  Z2,  Z3,  IV,  I  POE,  STRLT  > 

i  r  NON  .  'I  JT/IH.  HI,  fil-IC,  WE,  WN,  WT,  THETR,  PI,  RHOOO,  T,  CD,  Rfc,  VRHO,  VWX,  VWV 
'-•I HR.  VF,  VT,  INTfiLT,  ITIM,  CFRF.Q,  TNDl,  TND? 

I  >  "F  TH,  ITIM,  RND  INTFILT  FlS  INTERGER*4  VRR 
IH  1  r-2  IN'IRP,  IV,  I  PGP,  INCRFM 

INCFC-M  260 

.  F*-  I  WTirPP  FO  2>  I NCRF M~180G 
i  F  •  H  T  GT  ST  HI  T  >  GOTO  404 
IF  1 N  TRI  .7  .)  700,700,761 
MOOT  I H.  INCRFM 
INTRI  T  INTRLT*INCRF.M 
GO  TO  4G-! 

IF  ■  IN  i'RI  T  - 1 H  >  464,  702,  702 
h:  TRI  T  1 1,  f  i  "  T 

EE  I  Rl. T ;  <  fi! : T Rl . T -  SRI  >  /  <  H I  — SRI  > 

HE  1 1  SR3+LEI  Rl.T*<WF.-SR3> 

WUIT  SR-1+DEI  Rl  T*<WN-Sfi4> 

I'T  I T  SRG+DFLRLT  *  <  WT-SRG  ) 

CF'I  L  wrngi  cttt,  hfit,  wnit> 

IHF  r I T  TTl 

7.1  .  SR7+DELRI  T *<72-SR7> 

73  1  SR3+DEI  RLT *  <.  Z3-SR8  > 

*.t  r  r  srco+drlrlt* <  cfreg-sr^o > 

TNI  0  1 T  SR 2 1 +DFL RLT* < TNDl -SR71 > 

NIL-21  T  SR22+DFLRL  T*-' TND?-SR?2> 

V-  HO  I T  SR17+DRLRLT  +■  <  VRH0-SR1?) 

VWX  1  .  SR.1S+DEL  RLT  *  <  VWX— SR1S  > 

VNVI  T:  SR1.3+DFI  .RLT *<  VWV— SRJ  SO 
VT  I  i  SR10+DFI.RI  T*<V7-SR10> 

VF  IT  SR9+DELRLT* L VP-SR9 > 

P IMF: IT  SRI  i*<  03  01*RI/Sfill>**DELRL  T> 

R- 10 :  T  1  COG  0*  SRI  3*  <  < RH000/  SR1 3  )  **DEl.Rl.T  > 

I  7  I T  r  SR  12+r-Fl  R!  .7  *  T— SR12 1* 

CD  I  7  SHfi+r-'  L  RLT*  <  CD— SR14  > 

I RE  1 1-  SH 10+DFLRI  T*< RF-SR15 > 

RMC I T  SRI 6+DFLRLT*<  RMC-SR16  > 
fit  T- IN  PHI  T 
Tt  up-  ith 

IFCIPGF  LT  Gb>  GO  TO  33 
IPGF  3 

HR  I  TF.<  IV,  40:1  ) 

FuRMRTCIHl. 90X,  27H***R  M  S  NOISE  EkROR  IN***/3X,  119HTIME  RLT  . 


51 

52 

53 

54 

55 

56 

57 

58 

59 
68 
61 
62 

63 

64 

65 

66 

67 

68 
69 
78 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81. 
82 

83 

84 

85 

86 

87 

88 
89 


*  EWIND  NWIND  SPEED  DIR  PRFSS  T  DENSITY  Z  VEl.  Z  RCC  TDE 

*N  CF  DENS  PRES  TEMP  EW  NW/3X,  12SHZUL  U  METERS  M/S 

*  M/S  M/S  DEG  MB  K  GR/M3  M/S  M/S?  /CC 

*  /S  ♦♦♦♦♦•PERCENT*****  M/S  M/S  > 

33  IPGE:-IPGE+1 

WRITE  <4,  900)  fll  T,  TEMP,  RHOIT,  PIMBIT,  WEIT,  WNIT 
960  F0RM8T <  '  ' ,  6E10.  4) 

WRITE<  IV,  327>  INTfll  T,  WE  IT,  WNIT,  WTIT,  IHETIT,  PIMBIT,  ITIT,  RHOIT,  Z2IT, 
*Z3 IT,  TND2IT,  CFRQIT,  VRHOIT,  VPJT,  VTIT,  VWXIT,  VWVIT 
327  FORMRTOX,  17,  2F6.  1,  F7.  2,  2<I4,  F9.  5>,  F8.  2,  F6.  2,  2<1X,  E9.  4), 

*3F6.  1,  2F5.  1  > 

I NTRl  T- 1 NT8LT- 1 NCREM 
GO  TO  7©1 
404  S81-HI 
SR3-WE 

584-  WN 

585- WT 
S86=  THET8 
S8V-  Z2 
SR8--Z3 
SR9-VP 
SR10-VT 
SR11-PIMB 
SR12"T 
SR13-RH000 
SR14-CD 
SR15~RE 
SR16RMC 
SR17-VRH0 
SR18~VWX 
SR19-VWV 
SR20:::CFREQ 
SR21-TND1 
SR22--TND2 
RETURN 

tRSSM 

LIST 

*FORT 

END 
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tH  c -j  m 


*  $RSSM 
SCRAT 

I  INEAR  PROG  COMPUTF  L  INEAR  COEFFS  3©  NOV  73  R01  S  RFM3 :  LINEAR.  FOR 
‘1  CROSS 

0  NORX3 

6  NLS'I  C 

7  NLIST 

8  TFORT 


9 
,i  0 
:1  i 
;1  2 
1  3 
:1  4 
1  L. 
:i  6 
,:i : 
,3  C' 
19 
70 
PI 


SUBROUT  I NE  L.  I  NEAR  <  NF’TS  > 

COMMON  /COF-F/PH VI  <  51  ) ,  PXV2 <  35  > ,  P72  <  21  > ,  PZ1  <  51  > 
INTFGER+-2  NPTS..  M>  Kl,  K2,  I,  Ml 
7  p.-.  NPTS 

D-  <  <  TP-1 .  > *■  < TP  > *  < TP+1.  >  > /1 2. 

M:  CTP-l.  >/2. 
mi-  M+i 
P71  <  Ml.  > ~  ©. 

DO  26  I  1,  M 
Kl  Ml  I 
Li  M1.+  I 
P71  <K2>:  I/D 
20  P71 <K1  '  -°71<K2> 

RETURN 

TASSM 

LIST 

TKORT 

END 
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♦ASSM 


1 


2 

SCRfiT 

3 

RFAVG 

PROG  HftLF  SFC  AVERAGING  14  N0V149  R01  S  RFM3 :  RFAVG.  FOR 

4 

CROSS 

5 

N0RX3 

6 

NLSTC 

7 

NLIST 

8 

♦FORT 

9 

SUBROUTINE  RF.flVGdCNT,  IFLAG31 

10 

COMMON  T I ME  < 1 OO  > ,  X<100>,  V<100>,  Z<100>,  TMID<50>,  XMID<50>,  VMIDCSOl 

1.1 

ft  ZMID<50>,  XVM<50>,  VVM<50>,  ZVM<50>,  IX,  J,  N,  MIDMAX,  KMftXl,  KMAX2, 

12 

B  Cl,  NI  SI ,  NXV1,  NZ1,  ZXM<5S> 

13 

I NTLGFR*2  IX,  J,  N,  MIDMflX,  KMftXl,  KMAX2,  NIST,  NXV1,  NZ1 

14 

INTEGFR*2  I,  IFLAG3,  ICNT 

15 

IFLAG3— 1 

16 

ST=0. 

1.7 

SX=8. 

1.8 

SV=0. 

19 

SZ-0. 

28 

DO  1  1=1,5 

21. 

C 

DATA  FROM  THF.  BENCHMARK  RUN  TAPE  V388 

C 

AND  PROGRAM  8150SERNfl*BNCHMRK.  PMR 

23 

C 

PRODUCES  DATA  LABELED  "TIMF.  X  V  2". 

24 

C 

THF  FINAL  EDITED  DATA  CORRESPONDING  TO 

25 

C 

TO  "X"  <S1!>  IS  USED  FOR  CALCULATING 

26 

C 

VCJ),  AND  DAI  A  CORRESPONDING  TO  “V" 

27 

C 

<S2)  FOR  XCJ>. 

28 

C 

THIS  PROGRAM  WILL  FOLLOW  THIS  CONVENTION 

29 

C 

IN  THE  STATEMENTS  BETWEEN  LABELS  5  AND  1. 

38 

C 

DATA  READ  Sl  -X,  SX=V,  S2=Z 

31 

READ  <1,  5,  END=2>  T,  SI,  S2,  S3 

32 

5 

FORMAT  C4E1.S.  7) 

33 

ST=ST-*-T 

34 

SX=SX+S2 

35 

SV=SV+S1 

36 

1 

SZ=SZ+S3 

37 

TIME< J>=0.  2*ST 

38 

X<J>~  06096*SX 

39 

V<J>=.  06096*SV 

48 

2<J>=.  06096*52 

41 

RETURN 

42 

2 

IFL.AG3=0 

43 

RETURN 

44 

♦ASSM 

45 

LIST 

46 

♦FORT 

47 

END 
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*RS9H 


1 


4 


9 
j  0 
It 

13 

14 
i  ■: 
1 6 
1? 
.19 

19 
29 
2:1 
22 

24 

25 

26 

27 

28 
29 

20 


24 

25 

26 

27 

28 
29 


SCRfiT 

RFAVGT  PROG  THEORETICAL  AVGJNG  14  NOV  79  R03.  S  RKM3 : REAVGT.  FOR 
CROSS 
NORX3 
NLSTC 
NL  1ST 

IF-  ORT 

SUBROUT I NE  RFRVGT < I FLRG4 > 

COMMON  /THFO/ZT < 3 > ,  R2< 3  > ,  D2,  XT  <  3 > ,  RX <  3  > ,  DX,  VT  <  3 > ,  RV <  3  > , 

X  DV,  WX,  WV 

COMMON  T I ME < 3  80  > ,  X < 188  > ,  V  < 1 88  > ,  Z  < 180  > ,  TM I D  <  58  > ,  XM I D <  58  > ,  VM I D  <  58  > 
A  Zl'l IDF 08 > ,  XVM<58>,  WM(50)-  2VMF58),  IX..  J,  N,  MIDMRX,  KMRX1,  KMRX2, 

B  Cl,  NIST,  NX VI,  NZ1,  ZXMF50> 

INTEGERS  IX,  J,  N,  MIDMRX,  KMRX1,  KMRX2,  NIST,  NX VI,  NZ1 

INTEGERS  I,  IFL.RG4 

IFL.RG.p— 1 

ST  8.  0 

SX-8.  8 

SV-O.  8 

S2--0.  8 

DO  1  I- 1,5 

CALL  THFOT  <T> 

IF  FT.  EG!.  -999.  9>  GO  TO  18 

ST'-ST+T 

SX“ SX+XT <1 ) 

SVr- SV+VT  <  1  > 

1  SZ-SZ+ZT  FI  > 

TIMFF  J>=-0.  2+ ST 
X  <  J  )  ~8.  2+SX 
VFJV  0,  2+SV 
Z T )  ;"-G.  2*S2 
RETURN 
18  IFL.RG4-0 
RET  URN 

tRSSM 

LIST 

$FORT 

END 
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1 

*HSSM 

p 

SCRAT 

3 

91.  IDE 

PROG  MOVE  DATA  1  PT  14  NOV  79  R03 

4 

CROSS 

Li 

NGRX3 

6 

NLSTC 

7 

Nl.IST 

8 

*FORT 

3 

SUBROUTINE  SI  IDE <  IP,  ICNT,  IFl  RG5> 

10 

COMMON  TIMEC10G),  X<100>,  VC100>,  ?C10O>, 

13 

R 

1  ZMIDC50},  XVM < 50 > ,  VVM<50>,  ZVM<50>,  IX,  J 

•i ;  •■ 

B 

:  Cl,  NIST,  NXVl,  N21 , 2XM<50> 

3  3 

INTFGER*2  IX,  J,  N,  MIDMRX,  KMRX1,  KMRX2,  N3. 

14 

I NTEGFR+2  I,  IR,  IC,  ICNT,  IFLRG3,  I  FI  AGS 

15 

I  FLAGS-  -1 

3  6 

Ifl^KMRXl-IR 

17 

DO  3.18  I«l,  IR 

3  8 

IC-I+IB 

19 

XCI>;  X<  IC> 

20 

vci)=vac> 

21 

zc  i  >-z<  ic> 

22 

31© 

T 1  ME  C  I  '>  :-T  I MF  <  I C  > 

23 

IC-IR+l 

24 

DO  111  J-IC, KMRX1 

25 

13  3 

CALL  RE fl VG  < I CNT ,  I FLRG3  > 

26 

IF  Cl  FI  AGIO  114,  112,  113 

2  7 

3.14 

RETURN 

28 

3  32 

I FLAGS- 0 

29 

RETURN 

80 

3,13 

IFl  RGG--1 

31 

RETURN 

■?'£> 

1RSSM 

3  3 

LIST 

34 

tFORT 

35 

END 

RFM3 : SLIDF.  FOR 
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*R:-.sm 


SCRRT 

SI.IDET  PROG  MOVE  TDRTR  1  PT  14  NOV  79  ROl  S  REMS  :  SI  IDE T.  FOR 
CROSS 
N0RX3 

ni  etc 

NI..IST 

f-FORT 

SUBROUTINE  SL  IDFTCIR,  I  FI  RG6> 

COMMON  T I MF  <  1 OO  > ,  X  <  1 00  >  ,  V  <  100  >,  Z< 100 > ,  TM I D <  50  > ,  XII I D <  50  > ,  VM I D  <  50  >  , 
ft  ZM I D  <  5G  )  .•  XVM<59>,  VVM < 50  >  ,  ZVM<50>,  IX,  J,  1-4,  MIDMRX,  KMRX1,  KMRX2, 
e  Cl ,  NI  S  I ,  NXV1 ,  NZ1,  ZXI-1  < 58  > 

INTEGERS  IX,  J,  N,  MIDMRX,  KMRX1,  KMRX2,  NISI ,  NX VI,  NZ1 
INTF.GFR+2  I,  Ifl,  IC,  IFL.RG4,  I FLAGS 
IFL  RGS--  -1 
I  A:  KMRX1-IB 
DO  :'i.  'I  G  I  '•  j  ,  I A 
I  C:~  I  + 1 P 
X<  T  >•  X<IC> 
v ■  i V(IC) 

z< i>~z<ic> 

}  :i  0  TIMEt  I  ')■  TIMhCIC) 

I C-  T  R+l 

DO  111  J"  IC,  KMRX1 

111  CRl  L  RFRVGT <  IFL  AG4  > 

I F  <  I  FI  RG4  >  1 13,  112,113 
11  3  Rc.i  URN 

11 2  I FLAGS  0 
RFTURN 

*AS5M 

I  1ST 

•JfPORT 

END 
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CONSTR  PRI4T  Pt.OT  14  NOV  79  R01 


4  FORT 


SCRAT 
PROG 
CROSS 
NCR's  3 
NLS1 C 
NLIST 


SUBROUT  INF  TARCFSAL.T,  IV> 
INTEGER*2  I,  IV,  J 
D I  MENS  I  ON  DEN  <  5,  8  ) ,  W I MD  <  5,  8  ) 
DATA  ZERO/’  —  -/ 


RFM3 :  TAB.  FOR 


•13 

DATA  DEN 

/  3*4  HO.  01 

9 

4H8. 

47 

9 

4H0.  94, 

14 

A 

2*4H0.  01 

9 

4H0. 

02 

9 

4H0.  51, 

4H0.  95 

) 

15 

X 

2*4H0.  81 

9 

4H0. 

03 

9 

4H0.  59, 

4H0.  96 

) 

16 

X 

4H0.  01 

* 

4H8.  ©4 

9 

4H0.  47, 

4H0.  72 

9 

4H0. 

98 

) 

17 

X 

•4H0.  01 

9 

4H0. 

22 

.. 

4H0.  84, 

4H0.  97 

9 

4H0. 

99 

9 

18 

V 

A 

4H0.  22 

9 

4H0. 

52 

) 

4H0.  90, 

4H0.  98 

9 

4H0 

99 

9 

19 

yt 

4H0.  42 

9 

4H8. 

80 

9 

4H0.  96, 

4H0.  99 

9 

4H0. 

99 

20 

X 

4H0.  75 

9 

4H0. 

90 

) 

4H0.  99, 

4H0.  99 

9 

4H0. 

99 

/ 

21 

DATA  WIND 

/  4H0.  81 

9 

4H0. 

01 

9 

4H0.  62, 

4H8.  14 

9 

4H0. 

58 

9 

22 

X 

4H0.  01 

9 

4H0. 

01 

) 

4H0.  03, 

4H0.  17 

9 

4H0. 

63 

9 

23 

X 

4H0.  01 

9 

4H6. 

01 

) 

4HG.  06, 

4H0.  24 

9 

4H0. 

66 

9 

24 

.X 

4H8.  01 

.. 

4H0. 

03 

) 

4H0.  18, 

4H0.  85 

9 

4H0. 

90 

9 

25 

X 

4H0.  01 

.. 

4HG. 

16 

.. 

4H0.  78, 

4H0.  97 

9 

4H0. 

98 

, 

26 

X 

4H0.  02 

9 

4H0. 

40 

4H0.  98, 

4H0.  99 

9 

4H0. 

99 

9 

27 

X 

4H0.  57 

9 

4H0. 

76 

9 

4H0.  99, 

4H0.  99 

9 

4H0. 

99 

9 

28 

v 

A 

4H0.  96 

9 

4H0. 

99 

9 

4H6.  99, 

4H0.  99 

9 

4H0. 

99 

/ 

100  FORMAT < 1H1,  20X,  68HRAT 10  OF  AMPL  ITUDE  OF  SMOOTHED  DENSITY  WAVE  TO  A 
AMPLITUDE  OF  ORIGINAL ./32X, 46HWAVE  AS  A  FUNCTION  OF  ALTITUDE  AND  WAV 
2EI  ENGTH.  /29X,  58HC ALTITUDF  AND  WAVFL  ENGTH  MEASURED  IN  KILOMETERS.  ) 

l//> 

101  FORMAT C20X,  62HRATI0  OF  AMPl ITUDE  OF  SMOOTHED  SINUSOIDAL  WIND  TO  AM 
1PLITUDE  0F/23X,  55H0RIGINAI  WIND  AS  A  FUNCTION  OF  AL  TITUDE  AND  WAVE 
2LENGTH.  /29X,  5©H<  AIT  ITUDE  AND  WAVELENGTH  MEASURED  IN  KILOMETERS.  ) 

3//  ) 

a  02  FORMAT C46X,  SHALT I TUDE//26X,  57HX  180  X  90  X  80  X  70  X  60  X 

1  50  X  40  X  30  X/26X,  57<1HX)/20X,  7HW  X,  8C6X,  1HX)/20X,  7HA. 

2  1  X,  SCA6,  1HX>/20X,  7HV  X,  8C6X,  1HX)/20X,  7HE  2  X,  8C R6,  lHX)/20 

3X,  7HL  X,  8C6X,  1HXV20X,  7HE  5  X, 8CA6,  1HX)/20X,  7HN  X,  8C6X,  1 

4HX)/20X,  7HG  10  X,  8CA6,  1HX)/26X,  7HT  X,  8C6X,  1HX)/20X,  7HH  20  X, 

08 C  A6,  1HX > ////////// > 

WR I TE  <  I V,  1 00  ) 

IFCFSAl.T.  GT.  70.  )  GO  TO  20 

10  WR  I  TEC  IV,  10?)  CZERO,  ZERO,  ZERO,  ZERO,  CDENC  I,  J),  J=5,  8),  1=1,  5) 
WRITECIV,  181) 

WR  I  TE  C  I V,  1 02  )  <  ZERO,  ZERO,  ZERO,  ZERO,  C  W I ND  C  I ,  J  ) ,  J-  5,  8  ) ,  I  - 1,  5  ) 

RETURN 

20  IFCFSAL  T.  GT.  90.  )  GO  TO  30 

WR  IT  E  C  I V,  102  )  <  ZERO,  ZERO,  ZERO,  C  DEN  <  I ,  J  ) ,  J"4,  8  ) ,  1=1,5) 


WRITE  CIV/ 101 > 

WRITF.C  IV,  102)  CZERO,  ZERO,  ZERO,  CWINDC  I,  J>,  4,  8>,  1^1,  5> 

RETURN 

IFCF.SfiLT.  GT.  100.  >  GO  TO  40 

WR I  Te:  C  I V,  102  >  C  ZERO,  2ER0,  <  DFN  C  I ,  J  > ,  J--3,  8  > ,  I  -1,  5  > 

WRITE CIV,  101> 

WRITECIV,  10?>  CZERO,  ZERO,  CWINDC I,  J>,  J-2,  8>,  1-1,  5> 
RFTURN 

IFCESRl.T.  GT.  110.  >  GO  TO  50 

WR  I  TE  C I V,  1G2  >  <  ZERO,  <  PFNC  I ,  J  > ,  J=  2,  8  > ,  I  - 1,  5  > 

WRITE*  CIV,  101  > 

WR  I  TE  C  IV,  1 02 !)  < ZERO,  CWINDC  I,  J),  J--2 ,  8  > ,  I  '  l,  5> 

RETURN 

WR  I TF  <  I V,  1 02  >  C  C  DEN  <  I ,  J  > ,  J- 1 ,  8  > ,  I  1,  5  > 

WRITECIV,  10i> 

WR  I  TE  C  I V,  1 02  >  C  C bi I  NL> C  I,  J>,  J-l,  8),  I ~1,  5 > 

RETURN 


TH'-.SM 


*FORT 


$RSSM 


1 

2 

3 

4 


8 

9 

10 

1:1 

12 

13 

14 
1C. 
16 

17 

18 
.19 
28 
21 
22 

23 

24 
2’,. 
26 

27 

28 
29 
38 

31 

32 

.33 

34 

3C 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


SCRRT 

THfc'OT  PROG  COMP  THEO  Tkfl JKC  14  NOV  79  R01  S  RFM2 : THEOT.  FOR 
CROSS 
N0RX3 
NLSTC 
Nl.IST 

i-FORT 

SUBROUT INF  THFOTCT) 

INTEGERS  II,  K,  U,  IU 

COMMON  /CONST /  RG,  PHI,  HMSL,  RMS,  ZB,  DI,  VB,  RM,  GSRG,  RR 

COMMON  /THEO/  2<3>,  RZ<3>,  DZ,  X<2>,  RX<3>,  DX,  V<3>,  RV<3>,  DV,  WX,  WV 

DIMENSION  fl<4>,  B<4>,  PZ<2>,  PX<2>,  PV<2>,  QZ<2>,  QX<2>,  QV<2> 

DRTR  R,  B/0.  0,  O.  5,  &.  5,  1.  0,  1.  0,  2.  0,  2.  0,  1.  0/,  IU/0/ 

IF<  IU.  NF..  0>  GO  TO  2 

T~0. 

iu=-i 

2  DO  1  K--1,  2 
PZCK>^0.  0 
PXOO--0.  0 
PV<K>"0.  0 
QZCK>  =  0.  0 
QX<K>-0.  0 

1  GV<K)r-0.  0 

H— 0.  1 
DO  6  U-  l,  4 

DO  7  K-  1, 2 

RZ  C  K.  >  ■-■2.  (K)+h<U  >  +PZ  <  K  > 

RX  <  K  >  ~X  <  K  >  +R  <  U  )  *PX  <  K  > 

7  R V C K > -  V < K >  +fl < U 1  *PV <K> 

DO  S  K— 1,  2 

CALL  DRVT <K',  Rl,  R2,  R3,  R4,  R5,  R6,  Il> 

PZ <K>«D?*H 
PX  <  K  >  =~DX*H 
PV<K>--DV*H 

Q2  <  K  >  "-QZ  <  K  >  +B  <  U  >  *P2  <  K  > 

QX <  K  > ~QX <  K > +B(U> *PX  <  K  > 

QV  <  K  > = •  Q  V  <  K  > +B  <  U  >  *P' V  <  K  > 

6  CONTINUE 
DO  8  K"l,  2 
Z<K>-~Z<k>+QZ<K>/6. 

X  <  K  > “X  <  K  >  +QX  <  K  >  /£. 

8  V<K>-V<K>+QV<K>/6. 

T  "T  +H 

TRt  Tr--2<1  >+<X<l >*X<1 1+VC1 >*V<1 >  > / < 2*RG >  +HMSL 
IFCTR!  T.  GT.  55800.  )  GO  TO  10 
T— 999.  9 
IU=0 

10  RETURN 
$RSSM 

LIST 

*FORT 

END 
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*RSSM 

SCRAT 

TIFALL  PROG  CHK  SPH1  COLLAPSE  14  NOV  79  R01  S  RF.M3 :  TIFALL.  FOR 
CROSS 
N0RX2 
Nl.Sl  C 
NL  1ST 

♦FORT 

SUBROUTINE  TIFftt.L  <HI,  TFT,  I  FT,  ZVM,  HI2,  COLRPS,  IPOE,  IV> 

INTEGERS  IFT,  IH2,  IPGE,  IV,  JFT 

COMMON  /TF/TFT1,  TFT2,  TFT3,  TFT4,  TFT5,  TFT6,  TFT?,  TFT8,  TFT9,  TFT0 
DIMENSION  ZVHC50> 

GO  TO  CZ99,  302,  306,  300,  464,  312,  325,  360,  356,  362,  351  >,  IFT 

302  IFCZVM<2>.  LT.  ZVMC1>>  GO  TO  303 

412  IFT-5 
■JFT=1 

WRITE CIV,  41 1> 

411  FORMAT v20X, 58HRP0GEF  NOT  KNOWN.  TIME  OF  FALL  TEST  NOT  USED  ABOVE  5 

+5  KM.  > 

IPGF-  IF'GE+1 
GO  TO  399 

303  IF  CHI -779O0.  >  350,304,304 

304  IF  CZ\'MCl>+200.  >  412,305,305 

305  IFT=3 

GO  TO  399 

306  IF  C ZVM <  1  > -ZVM C  2  >  >  412,412,307 

307  IF  C  ZVM  Cl  >+20O.  >  308,308,399 

308  IFT=5 

HI2  - HI 21000.  +2.  1 

IH2-HI2 

HI3-1H2 

IF  CHI2-HI3-.  5>  309,310,310 

309  HI2=IH2 
GO  TO  311 

310  H72IH2+1 
I  H2-H  1 2 

311  WRITE. CIV,  4170  IH2 

413  FORMAT C1H  ,  20X,  16HBALL.00N  APOGEE  =,  14,  4H  KM.  > 

IPGE— IF’GF+l 

IFT.  6 
RETURN 

312  IFCTFT2  EG.  0.  01RETURN 

IF  CH 12-100.  >  300,301,313 

313  IFCHI2— 125.  >  314,318,319 

314  IF  CHI 2-115.  >  315,316,31? 

380  IFT-7 

GO  TO  325 
301  T  TEST-29. 

GO  TO  323 

315  TTEST-25.  +4.  *<115.  -HI21215. 


51 

GO  TO  323 

52 

316 

TTEST-25. 

S3 

GO  TO  323 

54 

317 

TTFST=24.  +<125.  -HI2)X10. 

55 

GO  TO  323 

56 

3d  8 

TTFST-24. 

57 

GO  TO  323 

58 

319 

IF<HI2-140.  )  320,  321,  322 

55* 

320 

TTEST=22.  +3.  *<140.  -HI 2) 215. 

68 

GO  TO  323 

61 

321 

TTEST-22. 

62 

GO  TO  323 

6.7 

322 

TTEST-21  +<150.  -HI2)/10. 

64 

323 

IF < TFT.  I  T.  TYEST+3.  .  AND.  TFT.  GT.  TTEST-3.  > 

GO 

to  : 

324 

65 

WRITF< IV,  41 4 > 

66 

414 

FORMAT < 20X,  76HT IMF  OF  FAIL  TEST  BETWEEN 

80 

AND 

70 

KM. 

67 

*LLOGN  DID  NOT  INFLATE.  > 

68 

IF<COl.APS.  EQ.  0.  )  COLAPS-70000. 

69 

IPGE--IPGF+1 

70 

IFT-7 

71 

RETURN 

72 

324 

WRI1  F< IV,  415) 

77 

415 

FORMAT < 28X,  32HBAU.00N  STILL.  INFLATED  AT 

70 

KM. 

) 

74 

IPGE=IPGE+1 

75 

IFT-7 

76 

RETURN 

77 

325 

IF<TFT3.  EQ.  0.  0)  RETURN 

78 

IF < TFT.  LT.  52.  .  AND.  TFT.  GT.  44.  )  GO  TO  326 

79 

WRITF< IV,  416) 

88 

416 

FORMAT < 20X,  57HT IMF  OF  FALL  TEST  BETWEEN 

70 

AND 

60 

KM. 

81 

Al  APSE) 

82 

IF< COLAPS.  FQ.  0.  )  COLAPS=60000. 

83 

IPGE=IPGF+1 

84 

IFT-8 

85 

RETURN 

86 

326 

WRITE<  IV,  417) 

8? 

417 

FORMAT < 20X,  32HBAI  . LOON  STILL  INFLATED  AT 

60 

KM. 

) 

88 

IPGF-IPGE+1 

89 

IFTr-8 

98 

RETURN 

91 

350 

I  FT  =4 

92 

IF<ZVM<1).  LT.  -150.  )  GO  TO  355 

9.3 

IFT=11 

94 

RETURN 

95 

351 

IF<ZVM<2).  GT.  ZVM<1>)  GO  TO  412 

96 

IF<ZVM<1).  GT.  -150.  )  RETURN 

97 

IFT=8 

98 

H 1 2=H 1X1000.  +1.  5 

99 

IH2=HI2 

100 

HI3-IH2 

INDICATES  BA 


INDICATES  CO 
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3.01  IFCHI2-HI2-.  5>  352,353,353 

I©?:  332  H 1 2=  I H2 

103  GO  TO  354 

104  353  HI2”IH?+1 

103  IH2-HI2 

106  354  WRITE CIV,  43  3  >  IH2 

18?  I FGE- IPGE+1 

IT'S  GO  TO  360 

105*  355  IPC7VMC1X  L.T.  -210.  >  GO  TO  432 

118  TFT  9 

3.11  RETURN 

112  356  I F  C  ?\  ;'K‘  2 > .  GT.  ZVM C 1  >  >  GO  TO  412 

113  IF <  ZVMC1  >.  GT.  -210.  )  RETURN 

114  H 1 2=H  1/1000.  +2.  1 

1.15  IH2---HI2 

116  HI 3- I HR 

11?  IFCHI2-HI3-.  5>  35?>358,  358 

118  35?  HI2:-  IH2 

119  GO  TO  359 

120  338  HI2-IH2+1 

121  IH2-H12 

1 22  339  WR I  TEC  IV,  41 3 >  IH2 

123  IPGfc.— IPGF+1 

124  360  IFCTET4.  EO.  0.  0>  RETURN 

125  I  FT -10 

126  IF < TFT  LT.  39.  .  AND.  TFT.  GT.  31.  >  GO  TO  361 

12?  UR 7 TEC  IV,  41  8 > 

128  41.8  FORMAT C20X,  5?HTIME  OF  FALL.  TEST  BETWEEN  60  AND  55  KM.  INDICATES  CO 

129  ALAPSE  > 

1.28  I FC  COL. APS.  EQ.  0.  >  COL APS--55000. 

131  I PGT— IPGE+1 

132  RETURN 

1 33  361  WR I TF.  <  I V,  41  9  > 

134  419  FORMAT C 2©X,  32HBALL00N  STILI.  INFLATED  AT  53  KM.  > 

135  IPGF  --  IPGE+1 

136  RETURN 

13?  464  GO  TO  <362,  365,  36?,  369,  3?1  >,  JFT 

138  362  IFCTFT5.  EQ.  0.  0>  RETURN 

139  IFT-5 

140  JFT -2 

141  IF C TFT.  LT.  54.  .  AND.  TFT.  GT.  44.  >  GO  TO  363 

142  WRITE  CIV,  428 > 

143  420  FORMAT < 20X,  57HT IMF.  OF  FAI  L  TEST  BETWEEN  55  RNI>  50  KM.  INDICATES  CO 

144  ALAPSE  > 

145  IFCCOL APS.  EQ.  0.  >  COLAPS^50000. 

146  I PGE- IPGE+1 

14?  RETURN 

148  363  WRITE CIV,  421 > 

149  421  FORMAT  C  20X,  32HBAL.L00N  STILL  INFL.ATED  AT  50  KM.  > 

150  I PGF= IPGF+1 
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151 

15? 

153 

154 
150 
156 
15? 

158 

159 
16ft 
161 
162 

163 

1 64 

165 

166 
1.6? 
168 
169 
17ft 
171 

177 

173 

174 

175 

176 
17? 

178 

179 
18ft 
181 
18? 

183 

184 

185 

186 

187 

188 
189 
19ft 
1.91 
197 

193 

194 

195 

196 

197 

198 

199 
2ftft 
2ftl 
2ft? 
2ft3 


399  RETURN 

365  IF  <TFT7.  LF.  0.  6)  GO  TO  399 
JFT--3 

IPGE-IPGF+l 

IF < TFT.  LT.  74.  .  AND.  TFT.  GT.  60.  )  GO  TO  366 
WRI7E<IV,  422) 

422  FORMAT  <  20X,  57HT I  ME  OF  FAIL  TF.ST  BETWEFN  50  AND  45  KM.  INDICATES  CO 
♦LAPSE  > 

IF<CQL APS.  EG.  0.  0)  COLAPS=45000. 

GO  TO  399 

366  WRI  TF/IV,  423) 

423  FORMAT < 20X,  32HBALL00N  STILi  INFl.ATFD  AT  45  KM.  ) 

GO  TO  399 

36?  IF  <TFT8.  LF.  0.  0)  GO  TO  399 
JFT  ~-4 

IPGE-IPGF+1 

IF<TFT.  LT.  106.  .  AND.  TFT.  GT.  84,  )  GO  TO  368 
WRITE CIV,  424) 

424  FORMAT < 20X,  57HTI ME  OF  FALL.  TEST  BETWEEN  45  AND  4ft  KM.  INDICATES  CO 
♦LAPSE  ) 

IFCCOLAPS.  EG.  0.  0)  COL. APS- 40060. 

GO  TO  399 

368  WRITE/.  IV,  425) 

425  FORMAT <20X,  32HBAI..L00N  STILL  INFLATED  AT  40  KM.  ) 

GO  TO  399 

369  IF  <TFT9.  LE.  0.  6)  GO  TO  399 
JF7  ■■5 

IPGE-=IPGE+1 

IF  <  TFT.  LT.  156  .  AND.  TFT.  GT.  124.  )  GO  TO  370 
WRITE/  IV,  426) 

426  FORMAT <20X..  5/HTIME  OF  FAIL  TEST  BFTWEF.N  40  AND  35  KM.  INDICATES  CO 
♦LAPSE  ) 

IFCCOLAPS.  EG.  0.  0)  COLAPS-  35000. 

GO  TO  399 

370  WRJTE< IV*  427) 

427  FORMAT <20X,  32HBAI.L00N  STILL  INFLATED  AT  35  KM.  ) 

GO  TO  399 

371  IF  CTFT0.  LE.  0.  0)  GO  TO  399 
I  FT-- 1 

IPGE-IPGF+1 

IF < TFT.  LT.  232.  .  AND.  TFT.  GT.  198.  )  GO  TO  372 
WRITE/ IV,  428) 

428  FORMAT <20X,  57HTIME  OF  FALL.  TEST  BETWEEN  35  AND  30  KM.  INDICATES  CO 
♦LAPSE  ) 

IF < COLAPS.  EQ.  0.  0)  COL APS-30G00. 

GO  TO  399 
37?  WRITE<IV, 429) 

429  FORMAT < 20X,  32HBAL.L00N  STILL  INFLATED  AT  30  KM.  ) 

GO  TO  399 

♦ASSM 

LIST 

♦FORT 


END 
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#ASSM 


SCRFIT 

“• 

TJFAl  2 

PROG  CHK  SPH2  COLLAPSE  14  NOV  79  R01  S  REM3  :  T I F AL  2.  FOR 

4 

CROSS 

N0RX3 

6 

NLSTC 

l’’ 

Nl.IST 

9 

*F0RT 

9 

SUBROUT  INF  TIFAL  2<HI,  TFT,  IFT,  ZVM,  HI2,  COLAPS,  IPQF,  IV> 

J  0 

INTE6FP+2  IFT . JFT,  IH2,  IPGE,  IV 

1 

COMMON  .'TF/TKT1,  TRT2,  TFT3,  TFT4,  TKT5,  TRT6,  TFT7,  TFT8,  TFT9,  TFT0 

1  * 

D  I  MENS  I  ON  ZVM < 58) 

4  ' 

GO  TO  <399,  302,  306,  300,  464,  212,  320,  360,  306,  362,  351  >,  IFT 

14 

202 

I F  <  ZVM  <  2 > .  LT  ZV M<J>>  00  TO  302 

10 

412 

IFT=5 

1  6 

JFT  1 

1? 

WRITfc  <IV,  411) 

1.8 

411 

FORMAT <P6X,  08HAP0GEE  NOT  KNOWN.  TIME  OF  FALL  TFS1  NOT 

USED  ABOVE 

19 

* 

0  KM  > 

■ti 

IF'GF"-  IPGF+1 

21 

GO  TO  399 

i-.-'v 

Ji-J 

IF  CHI  -779O0  )  350,204,304 

22 

204 

IF  < ZVM <  1  > +200.  >  412,300,200 

24 

20' . 

IFT:  3 

25 

GO  TO  299 

28 

3G6 

I F  <  ZVM <  1  > -ZVM  <  2  >  )  412,  412,  307 

2  V 

207 

IF  <ZVM<l>+2©0.  >  208,308,299 

28 

206 

IFT--5 

29 

HI  2  -HI/1000.  +2.  1 

30 

IH2  :HI2 

21 

HI 3- I HP 

22 

IF  <HI2-HI3-.  5>  309,310,31© 

2 2 

209 

HI  2=  JH? 

24 

GO  TO  311 

210 

HI  2--  -IH2+1 

26 

IH2- HI2 

27 

211 

WRITE<  IV,  41.3)  IH2 

28 

413 

FORMAT  <1H  ,  20X,  16HBALL00N  APOGEE  «=,  14,  4H  KM.  > 

29 

IFgF=IPGE+1 

48 

IFT -6 

41. 

RETURN 

42 

312 

IF<TFT2.  F.Q  0.  0)  RETURN 

42 

TTEST-23. 

44 

IF<TFT.  LT.  TTFST+O.  .  AND.  TFT.  GT.  TTF.ST-5.  >  GO  TO  324 

40 

WRITfc  < IV,  414) 

46 

414 

FORMAT<20X,  76HTIME  OF  FALL  TEST  BE TWEEN  80  AND  70  KM. 

INDICATES 

47 

♦LLGON  DID  NOT  INFLATE.  > 

48 

IF < COL  APS.  FQ.  0  )  COl.APS-70000. 

49 

IF'GE-IPGE+1 

oo 

IFT-"’ 

f 


51 

RETURN 

52 

124 

WRITE <  IV,  415) 

5.1 

415 

FORMAT < 20X,  32HBAL LOON  ST11.L  INFLATED 

AT 

70 

KM. 

54 

IPGE-IPGE+1 

55 

IFT=-7 

56 

RETURN 

5? 

125 

IF<TFT2.  E Q.  0.  0)  RFTURN 

58 

IF < TFT.  LT.  46.  .  AND.  TFT.  GT.  38.  )G0  TO  326 

59 

WRITEdV,  416) 

60 

416 

FORMAT<20X,  57HTIME  OF  FALL  TEST  BFTWEFN 

70 

AND 

61 

Al  APSE) 

62 

I F  <  COLAPS.  EQ.  0.  )  COLAPS=60OO0. 

61 

IPGF-IPGE+1 

64 

I  FT  --8 

65 

RETURN 

66 

126 

WRITE  <  IV,  417) 

6? 

417 

FORMAT < 20X,  32HBAL.L00N  STILL  INFLATED 

AT 

60 

KM. 

68 

IPGF-IPGE+1 

69 

I  FT -8 

70 

RETURN 

71 

ISO 

I  FT  "4 

77 

IFCZVMCl).  LT.  -150.  )  GO  TO  255 

71 

IFT=11 

74 

RETURN 

75 

151 

IFC7VMC2).  GT.  ZVM<1))  GO  TO  412 

76 

IFCZVMCl).  GT.  -150.  )  RETURN 

77 

IF  T=8 

78 

H 1 2-H  I  /1 000.  +1.  5 

79 

IH2=HI? 

80 

HI3-IH2 

81 

IF<HI2— HI3-.  5)  352,  353,  352 

87 

152 

HI2-IH? 

81< 

GO  TO  354 

84 

353 

HI7~IH2+1 

85 

I H2-H 1 2 

86 

154 

WRITE<IV,  413)  IH2 

87 

IPGF-IPGE+1 

88 

GO  TO  360 

89 

355 

I F < Z VM <  1 ) .  LT.  -210.  )  GO  TO  412 

90 

IFT=9 

93 

RFTURN 

92 

356 

IF<ZVM<2).  GT.  ZVM<1)>  GO  TO  41.2 

91 

IF<ZVM<1)  GT.  -210.  )  RETURN 

94 

H 1 2=H  1 71000.  +2.  1 

95 

IH2-HI2 

96 

HI3=IH2 

97 

IFCHI2-HI3-  5)  357,  358,  358 

98 

357 

HI2«IH2 

99 

GO  TO  359 

100 

258 

HI?«IH2+1 

INDICATES  CO 


129 


a  oi 
3  ox 
io< 
104 
10‘  • 
1  Or: 
10? 
It. 
109 
13  0 
13  1 


IH2=HT£ 

359  WR I TF.  <  I V,  43  2)  I H2 
IF‘GF:-  IPGE+1 

360  IFCTFT4,  FQ.  0.  0>  RETURN 
I FT “10 

IF ''TFT  L.T.  35.  .  FIND.  TFT.  GT.  2?.  >00  TO  363 
WRITE CIV,  418 > 

43  8  FORMAT C2GX,  57HTIME  OF  FALL  TEST  BFTWEEN  60  AND  55  KM. 
HL  AF’SF  > 

IF  COOL  APS.  EG)  0.  >  COLAPS-=550G0. 

IPGF“ IPGE+i 


13.  2 

RETURN 

1 3  3 

363 

WRITE1'  IV,  41 9 > 

1.1  4 

43  9 

FORMAT C28X, 32HBAI  LOON  Si  ILL  INFLATED  AT 

115 

I  ?GF.  ••  I PGF+1 

116 

RETURN 

3.1? 

464 

GO  TO  <362,  365,  36?,  369,  373.  >,  JFT 

118 

H>6<? 

TFCTH5.  EG!  0.  0>  RETURN 

1:1  9 

i  FT  : 

1 20 

JFT- 2 

1.  Vi 

IF  OFT  IT.  47.  .  AND.  TFT.  GT.  36.  >G0  TO  363 

122 

WR TEC  I V,  426) 

i-.‘  > 

426 

FORMAT <2©X,  57HTIME  OF  FAIL  TEST  BETWEEN 

3  24 

H1  hF'SE  > 

x,':  5 

I F  C  COl.APS.  FQ.  0.  >  COL  APS -50000. 

126 

I PGF.-  IF'GF+l 

12:? 

RETURN 

i  2'2 

36? 

WRITE'.:  IV,  421) 

129 

421 

FORMAT C20X,  32HBAL1.00N  STILL  INFLATED  AT 

1-0 

I PGF  IPGE+1 

1  J 

PE TURN 

1.  2 

36 

1“  cTFT,  L.E.  0.  0 >  GO  TO  399 

13 

.if  r  3 

3  <4 

I  PGF-  IFGF+1 

3  25 

iP  '"'‘FT  LT  65.  .  AND.  TFT  GT  50.  >  GO  TO  366 

136 

NR  I  TE  •:  IV,  422) 

.1  </ 

422 

FORMAT <20X.  57HTIMF.  OF  FAI  L  TEST  BETWEEN 

138 

♦I.AF'SE  > 

J  ■*-< 

IF  .COLAPS  EG!  0.  0>  COL. APS1  -45000. 

140 

GO  1  u  399 

1  4  J 

j  66 

WPI'ir  C  IV,  423 > 

i  42 

42? 

FORMAT <20X,  32HBAU  OOH  STIL  L  INELATFD  AT 

14  3 

GO  TO  399 

144 

~  6 

IF-  TFT8  LE  0.  0>  GO  TO  399 

145 

JFT  =4 

)  46. 

IPGFIPGE+l 

J  4? 

IF^TRT  LT  93.  .  AND.  TFT.  GT.  72.  >  GO  TO  368 

WRITE  CIV,  424  > 

J  44 

424 

FORMAT  <  20X,  57HTIME  OF  FAIL.  TEST  BETWEEN 

1"  0 

*1  APSE  > 

INDICATES  CO 


INDICATES  CO 


t 

i 


INDICATES  CO 


INDICATES  CO  * 
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151  IF < COt  ftPS  EG  0.  0)  COLfiPS~40000. 

152  GO  TO  399 

153  368  WRITE CIV,  425  > 

154  425  FORMflT<20X,  32HBftl.L00N  STILl.  INFLATED  FIT  40  KM.  > 

155  GO  TO  399 

156  369  CONTINUE 

157  371  CONTINUE 

158  GO  TO  399 

159  #RSSI1 

160  LIST 

161  *FGRT 

162  END 
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..  •  '.Sjt* 


*HSSM 


TP OB IN 


;  FO 


SCR  ft  T 

PROG  THEORECT  SRH  TRftJ  2G  NOV  79  R6J  S  REM2 : TROBJN.  FOR 
CROSS 
NORXS: 

Nl  ST C 
NLIST 

SUrsROU T I J-4F  TRORII'-KGS,  N2MID,  KSW2,  KSW,  IF'GF,  IV> 

If-'TF.GER+S  I  CD,  IFLftG,  I  FLftG-1 ,  IFLRGB,  I PGE .  ISN,  IV 
•  :  FGF.P+2  KSM,  KSW2,  I-4PRT,  N2MID,  JB,  I  MX,  JMX 

T.N'  -79ER+-2  IX,  -J,  N,  MIDMftX,  KMRX1,  KMRX7,  NIST,  NXV1,  N21 
'-,'TFGFP*?  I,  K,  NX  Vi?,  NZ2,  NX1I-1ID,  NX2MID,  NZiMID,  N22MID,  N2S1 
COMHO!-,'  .  'EXTkh/  NXV2,  NZ2,  NXiMID,  NX2MID,  1-471 H I D,  NZ2MID,  N2ST 
rr-;!'-HO'-4  /COFK7  PXVK51  >,  PXV205>,  PZ2<?1>,  P21  <51  > 

COMMON  /CONST/  RQ.  ftl.ft,  HMSL  ,  RMS,  ZB,  DIR,  VB,  RM,  GSRG,  RB 
■:  OMMON  /COR  1 2/  TOSL,  TOCL,  TOCL  SZ,  T0CLC7 
1  -i  '!».•;-<  /CON?/  ft  I  IK,  ft!  F'Hft,  PU,  RftD 
COMMON  /TftPL/  I  MX,  E!  I  ftSH  I  <  666  ) ,  B I  ftSMX  <  666  !) ,  B I  ftSI-JV '!  COO  !) , 

B I  ft  7R  v  BOO  > ,  E:  I  ftsT  *-.  SwO > ,  JMX,  B I RSMC  < S06  > 

COMMON  ,-■/  T I  ME  <  1 OO  > ,  X<1©6>,  V<106>,  7<10G>,  TMID<50>,  XMID<50>, 

Vi  i  I D < 5G  > ,  ZM I D < 56  > ,  XVM < 56 > ,  VVM  < 50  > ,  ZVM <  56  > ,  I X,  J,  N, 
MIDMftX,  KI-1RX1,  KMftX?,  Cl  -  NIST,  l-JXVl , 

:  1-421.  ZXM  <  56  > 

Tft  6 


,i  c .  '  1 
HZ’!  >  6 

!  Ir'RT  ■■  .1 

001  ->ii  52  J-  1,  KMftXl 
j  -  Hi  !....  ivErt'v'G  T I FL  .ftC’-l  ■' 

IF-'  JFL  RG4>6,  90,  € 

R  GO  TO  <?,  9>,  KSM2 
?  DO  8  K--1,  Ki-:ftX2 
r-4  K 

CROC  Pi  TON 

Cfti  .L  SI  .IDE  I  <2,  IFL  RGB  I1 
; F •"  IFL. RGB >8,  90,  8 

8  L.  Ui-47  INUE 
GO  TO  11 

9  DO  1 6  K-  1 ,  KMRX2 
N  -K 

CRLl  FIT0N2 
Cfil.L  SL. IDE  i  <2,  IFL  RGB!) 
IF1-.  IFLfiUb  )10,  96,  1W 
16  CONTINUE 

1..J.  XZ;  - 0.  0 

z:<x-  6.  0 

VC-  0.  6 
Z2;  6.  0 
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51 

GO  TO  <12,  15  >,  KSN 

%.*  f* 

12 

DO  12  1=1,  NXV2 

53 

22X  -  Z2X+PXV2  <  1  >*ZXM<  I  > 

54 

X2=X3+PXV2<  1 >+XVW<  I  > 

5!;. 

12 

V2"'V3+PXV2  <  I  >  *WM  <  I  > 

56 

NEND:  N2ST +N72  -1 

5? 

DO  14  I"N2ST,  NFND 

58 

K-I-N2S1+1 

59 

14 

Z3: :  73+P22  <  K  >  *ZVM  <  I  > 

66 

GO  TO  1 8 

61 

15 

DO  16  1=1, N22 

62 

16 

Z5:;  23+PZ2<  I  >*ZVM<  I  > 

67? 

WEND ' NXV2+M2ST  -1 

64 

DO  17  IN2ST,  NFND 

65 

K“ I -N2ST+1 

66 

Z3Xr-  Z2X+PX V2  <  K  > *7X14 <  I  > 

67 

X2=X3+PXV2 <  K  >  *XVM  <  I  > 

68 

17 

V3: ;  V3+PXV2  C  K  >  * VVM  <  I  > 

69 

18 

CONTINUE 

78 

X2;-X  VM  <  N2M 1 D  > 

7.1 

V2:- VVN  <  N2M I D  > 

1  .*J 

Z2X:-  ZXM  <  N2M I D  > 

7.“' 

2. 2  -ZVM  <  N2M I D  > 

74 

XI:  XMID<N2MID> 

75 

Vl:  VM  I D  <  N2U I D  > 

76 

Z1'ZMID<N2HID> 

7  7 

TIM--TMID<N2Mir>> 

78 

H I :  Z1  +  < X1+X1+V1 ♦VI  >/<2.  *RG >  +HMSL 

79 

OPHOR:  •  1 .  0+H I  /RG 

80 

UTHR- GS/ < RG*OPHOR*OPHOR+GPHOR > 

8:1 

IF?. 73.  LT.  —8.  0>  GO  TO  2050 

8? 

COR  X  ~  T  0  C  L  S  Z  •+=  7  7!  X + T  OS  L  *  V  2 

82 

CORV- T0CLCZ+Z2X-T0SL *X2 

84 

C0RZ::-T0CLSZ*X2-T0Cl.CZ+V2 

85 

GO  TO  ?19,  38  >,  ISM 

86 

19 

CRLL  flTNOS<HI,  T,  RHOO > 

87 

I  CD- 5 

88 

VT= ftBS<Z2>*SQRT<l.  6+<X3*X2+V3*V3>7 

y9 

PF:DI  fi+RHOO+VT *■  <  T+116.  4  >  /  <  O.  066601- 

9© 

fll1C=VT  +■  <  SORT  <  288.  16/T>>/340.  29265 

91 

CFtl.L  DRftGT  <  CD,  RMC,  RE,  I FLAG > 

92 

IF < I FLOG >34,  2666,  34 

97? 

24 

I  CD  1 

94 

CD1--CD 

95 

RHOGCi:  RHOO 

96 

RHO-RHOO 

97 

GO  TO  26 

98 

28 

DCD-CD 

99 

CD-2  6+CD-CDl 

.00 

CD1:  DCD 

P’l+RG)  )**2> 
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.X* 


n 


26  WTWO~  <:  :1 ,  G-VB+RHO/RMS  >  *WTHR 
BGIJVX:  -Xl+WTWu 

BGUVV  Vl+WTHO 
BGUVZ--  Zi+RG :>  +WTW0 

21  WDE  NO  IT  •  Z  S  X + C  0  R  2  -  El  0  U  V  Z 

22  WZ  WZ1 

HX-X2-  <  Z2X-MZ  >  *  <  X  3 + C  0  R  X  -  B  G  U  V  X  ) XWDF.NGM 
MV~  t‘2-  < Z2X-WZ > * C Y3+C0RV— BOUW  > XWDFNOM 
WZ1 :  < X1+WX+ V±*M Y  >  X  <  H I +RG > 

IF'  RF"~ *; HZ  L-WZ > .  GF,  O.  1>  GO  TO  22 
.■X  4:71 

VT :  SORT  <:  (. X2-WX >  **?+  <  V2-WV >*+2+  < Z2-WZ>**2 > 

F-Ho  Rl-i  ~  -•  <  -CORZ-23+WTHR*  <  RG+Z1  >  > /  <  CD+RR+VT  +  <  72-WZ  >  +WTHR*  <  RG+21  >  *VB  > 

GO  TO  23,  24 .■  ISM 

23  PI:  T+RU+RHO/RMk 
GO  TO  20 

24  P I  -  P I  T  +GF+RG**2+  <HI— HI  I  > *E XP< G.  5*RI.0G < RHO*RHOOO >  > / <  < RG+G.  5*  < H I +H 1 1 

H  '  >*+-2  > 

r  ~p  I  '+.R|,'ik/ RU*RHQ> 

20  Rl-iC  VT*<  SORT  <288.  1S/T>>X3*16.  29200 

kF  -0 I  R*RHO*VT *4  T +11  6.  4  >/<6.  600061458*SQRT < T **3 >  > 

COG ~ CD 

CRU  DRRGT  CD.*  RHC.-  RF.*  I  FI  RG.> 

T  t-  >;  |  p  |  R IS !'  j;-! 

32  ICD---J 

1  F  Nr'RT.  F i 01  UU  TO  1^6 
CO"CDG+  <  CO— CD©  > /3. 

I F  RRS  •:  RHO-RHOO > /RHO.  LT.  0.  663333 >  GO  TO  28 
IF'  UPRT  GF.  5>  GO  TO  27 
.  ‘ R  RHgr  -RHO 

i4F-i  T  MF'RT+i 
GO  T 0  20 

27  RH'.i'  RHU+RH001/2.  0 
C  0  •••  i '  0 + C  0  0  /  2 .  0 
NPR7  0 

GO  TO  20 
20  RHOOO - RHO 
NF'kT -  G 
RHUU- ;  RHO 
P i  T  PI 
HI  I  HI 

T F < 72.  IT  -8.  CO  GO  TO  2606 
ISW  2 

CRLI.  RTMOStHF  TSS..  RSS) 

JF:-JF:+1 

IF '  JR,  GT.  660)  GO  TO  2666 
BIRSHI  •:  jb)  --hi 
E:  I RSI1C  <■.  JB  )  "RMC 
BIRSMXc.  JB MX 
BIRSMVC  JB)--WV 
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*RSSM 


1. 

2  SCRflT 

3  WflNGI.  PROG  POLRR  WIND  COORD  14  NOV  79  R01  S  RFI12 :  WRNGI  FOR 

4  CROSS 

0  N0RX3 

6  NLSIC 

7  NLIST 

8  f FORT 

9  SUPROUT I NF  WftNGL  <  THFTR,  WF,  WN  > 

1 8  i NTEGFR+2  I ,  NFND,  N2M I D 


1.1 

RftD=  07.  2907795 

i 

IF  <WF  >1003..  1662..  1061 

i3 

1 681 

I F  ('  WN  >  1 011,  101:1 ,1010 

14 

1G82 

IFCWN>1014,  1614,  1010 

IS 

1.083 

I F  >:  WN  >  101 3,  161 2,  1G1  2 

If- 

1818 

TH-2.  141  5927+RTRN  <WE/WN> 

1  ? 

GO  TO  1018 

1  8 

1011 

TH  4.  7123898+87 ftN  <-WN/WF> 

19 

GO  TO  1618 

P& 

1012 

TH  1.  G7G7963+RTRN  <-WN/WF.> 

2:1. 

GO  TO  1618 

4.01.3 

TH"-  AT  AN  <  WF/WN  > 

>.• *£ 

GO  TO  101.8 

PA 

161  4 

TH  8.  8 

25 

1818 

THE  f8:;  TH  >  +  C  R8D  > 

26 

RETURN 

27* 

1RSSI-1 

2o 

LIST 

29 

stFuRT 

7c  0 

END 
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